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ABSTRACT 

A subset of blazars are powerful TeV emitters, dominating the extragalactic component of the very high energy 
gamma-ray universe (E > lOOGeV). These TeV gamma rays generate ultra-relativistic electron-positron pairs via 
pair production with the extragalactic background light. While it has generally been assumed tha t the kinetic en- 
ergy o f these pairs cascade to GeV gamma rays via inverse Compton scattering, we have argued in lBroderick et al.l 
(120121 Paper I in this series) that plasma beam instabilities are capable of dissipating the pairs' energy locally on 
timescales short in comparison to the inverse-Compton cooling time, heating the intergalactic medium (IGM) with 
a rate that is independent of density. This dramatically increases the entropy of the IGM after redshift z ~ 2, with 
a number of important implications for structure formation: (1) this suggests a scenario for the origin of the cool 
core (CC)/non-cool core (NCC) bimodality in galaxy clusters and groups. Early-forming galaxy groups are unaf- 
fected because they can efficiently radiate the additional entropy, developing a CC. However, late-forming groups 
do not have sufficient time to cool before the entropy is gravitationally reprocessed through successive mergers — 
counteracting cooling and potentially raising the core entropy further. This may result in a population of X-ray 
dim groups/clusters, consistent with X-ray stacking analyses of optically selected samples. Hence blazar heating 
works different than feedback by active galactic nuclei, which we show can balance radiative cooling but is unable 
to transform CC into NCC clusters on the buoyancy timescale due to the weak coupling between the mechanical 
energy to the cluster gas. (2) We predict a suppression of the Sunyaev-ZeFdovich (SZ) power spectrum template 
on angular scales smaller than 5' due to the globally reduced central pressure of groups and clusters forming after 
Z ~ 1. This allows for a larger rms amplitude of the density power spectrum, erg, and may reconcile SZ-inferred 
values with those by other cosmological probes even after allowing for a contribution due to patchy reionization. 
(3) Our redshift dependent entropy floor increases the characteristic halo mass below which dwarf galaxies cannot 
form by a factor of approximately 10 (50) at mean density (in voids) over that found in models that include pho- 
toionization alone. This prevents the formation of late-forming dwarf galaxies (z < 2) with masses ranging from 
10 10 to 10 11 M Q for redshifts z ~ 2 to 0, respectively. This may help resolve the "missing satellites problem" in 
the Milky Way of the low observed abundances of dwarf satellites compared to cold dark matter simulations and 
may bring the observed early star formation histories into agreement with galaxy formation models. At the same 
time, it explains the "void phenomenon" by suppressing the formation of galaxies within existing dwarf halos of 
masses < 3 x 10 10 Mq with a maximum circular velocity < 60 km s -1 for z < 2; hence reconciling the number of 
dwarfs in low-density regions in simulations and the paucity of those in observations. 

Subject headings: BL Lacertae objects: general - galaxies: clusters: general - galaxies: formation - galaxies: 
dwarf - gamma rays: general - intergalactic medium 



1. INTRODUCTION 

Extragalactic relativistic jets are powered by accreting super- 
massive black holes (or in general the engines of active galac- 
tic nuclei, AGNs) and are able to carry an enormous amount of 
power out to cosmological distances. Blazars are a subclass of 
AGNs where the jet opening angle of typically ~ 10° contains 
our line-of-sight, allowing us to detect the Doppler-boosted ra- 
diation. Blazars are the dominant extragalactic source class in 
the TeV sky with currently 36 known objects out of 46 extra- 
galactic sources (of the remaining 10, 4 are radio galaxies, 2 
are starb urst galaxies, and 4 are n ot yet identified; for a re- 
view, see iHinton & Hofma nn 2009)0 Recent observations by 
the Fermi Space Telescope and ground based imaging atmo- 

1 For an up-to-date list/visualization of the extragalactic TeV sky, see 
http://www.mppmu.mpg.de/~rwagner/sources/ or 
http://tevcat.uchicago.edu/ 



spheric Cherenkov telescopes (H.E.S.S., MAGIC, and VER- 
ITASj3 demonstrated that most of the electromagnetic power 
is emitted in the gamma-ray band. The most extreme blazars 
achieve energies of up to 10 TeV, giving rise to the class of 
high-energy peaked BL Lac objects (HBL) while the some- 
what less efficient accelerators in intermediate-energy peaked 
BL Lac objects (IBL) are, in some cases, also able to reach en- 
ergies beyond 100 GeV. The emission mechanism is thought to 
be inverse Compton scattering of ultra-relativistic electrons in 
the jet giving rise to power-law energy spectra that increase as a 
function of energy and peak at the maximum energy that the ac- 
celerating process of the radiating relativistic electrons is able 
to deliver. The universe is not transparent to very-high energy 

2 High Energy Stereoscopic System, Major Atmospheric Gamma Imaging 
Cerenkov Telescope, and Very Energetic Radiation Imaging Telescope Array 
System. 
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gamma-ray radiation (VHEGR; E > 100 GeV), i.e., a beam of 
these energetic photons will necessarily produce electron and 
positron pairs off of the extragalactic background light (EBL), 
with typical mean free paths of VHEGRs ranging from 30 Mpc 
to 1 Gpc depending upon gamma-ray energy and source red- 
shift. The pairs produced by VHEGR radiation have typical 
Lorentz factors of 10 5 - 10 7 . 

There are only two possible ultimate destinations for the ki- 
netic energy of the pairs: the energy can either be channeled 
into lower energy (GeV) gamma-ray radiation (to which the 
universe is transparent), or heat the ambient medium with a 
partitioning factor that depends on the relative rates of the pro- 
cesses. The first process has generally been assumed to dom- 
inate the manner in which these pairs lose energy, almost ex- 
clusively through inverse Compton scattering the cosmic mi- 
crowave background (CMB) and EBL on a typical mean free 
path of (10- 100) kpc today. When the up-scattered gamma- 
ray is itself a VHEGR the process repeats, creating a second 
generation of pairs and up-scattering additional photons. The 
result is an inverse Compton cascade depositing the energy of 
the original VHEGR in gamma rays with energies < 100 GeV. 

There are, however, problems with this picture. First, the 
expected inverse Compton bump has not been seen in the spec- 
tra of luminous blazars around lOGeV. This could imply the 
existence of intergalactic magnetic fields that deflect the pairs 
out of our lin e-of-sight, hence reduc i ng the inverse Comp- 
ton emission (Neronov & Vovkl |2O!0t iTavecchio et all 120 Id 



1 201 UlDermer et alj201 UlTavlor et alj201 U lDolag et alj201 1 



iTakahashi et alj|2012t IVovk et alj|2012l) . Typical lower limits 
for magnetic fields would then range from 10~ 19 G to 10~ 15 G, 
depending on the assumed duty cycle of blazars and are dom- 
inated by void regions, which dominate a typical line-of-sight. 
The values at the upper end are of astrophysical interest in 
the context of the formation of galactic fieldfl Second, the 
spectral shape of the unresolved extragalactic gamma-ray back- 
ground (EGRB), that has been measured by Fer mi, exhibits 
a ste ep power-law at energies below 100 GeV dAbdo et al.l 
120101) . If blazars contribute substantially to the EGRB, we 
would expect to see a flattening toward high energies, due 
to the inverse Compton cascades, in conflict with the Fermi 
EGRB. Traditionally, the EGRB is then used to constrain the 
evol ution of the luminos i ty density of VHEGR sources (see , 
e.g.. iNarumoto & TotanH 120061: iKneiske & Mannheim! 120081 : 
llnoue & Totanill2009t lVentersll2010l) . Generally, it has been 
found that these cannot have exhibited the dramatic rise in 
numbers by z ~ 1-2 seen in the quasar distribution. That is, 
the co-moving number of blazars must have remained essen- 
tially fixed, at odds with both the large-scale mass assembly 
history of the universe, e.g., star formation history, and with 
the luminosity history of similarly accreting systems, e.g., the 
quasar luminosity density. 

In our first companion paper of this series dBroderick et al.l 
|20T1 hereafter Paper I), we argued that instead of initiating an 
inverse Compton cascade, the pairs dissipate their kinetic en- 
ergy locally, heating the intergalactic medium (IGM). We iden- 
tified a process that operates on a timescale fast in comparison 
to inverse Compton cooling, dominating the latter for luminous 
TeV blazars, i.e., for HBL and IBL blazars with an equivalent 
isotropic luminosity of L > 10 42 ergs _1 above 100 GeV. De- 
spite its dilute nature, the VHEGR-generated beam of ultra- 

3 After adiabatic contraction and a handful of windings, nG field strengths 
can be produced from an intergalactic magnetic field of ~ 1(T 15 G. 



relativistic pairs, propagating through the IGM, is susceptible 
to plasma beam instabilities. While the commonly discussed 
Weibel and two-stream instabilities are strongly suppressed by 
finite beam temperatures, these are special cases of a general 
filamentary "oblique" instability which is far mor e virulent, 
and strongly insensitive to finite temperature effects dBret et al. 
200 1 l20Q5t iBrei 120091: UTret et alJl2010b iLemoine & Pelletier 
20101) . 

If these instabilities saturate at a rate that is comparable to 
their linear growth rate, the beam kinetic energy is directly 
transferred to electrostatic modes which rapidly dissipate lo- 
cally, heating the IGM. If this scenario or an analogous, sim- 
ilarly efficient mechanism operates in practice, it necessarily 
suppresses the inverse Compton cascades and naturally ex- 
plains the absence of an inverse Compton bump in the TeV 
blazar spectra without invoking an intergalactic magnetic field. 
At the same time, it allows for a redshift evolution of TeV 
blazars that is identical to that of quasars without overproduc- 
ing the EGRB. In fact, for plausible parameters of TeV blazar 
spectra, it is possible to explain the high-energy part of the 
EGRB (Paper I) without the need to appeal to exotic ph enom- 
ena (e.g., dark matter annihilation, iCavadini et alj|201 11) . 

By dissipating the pairs' energy into the IGM, plasma in- 
stabilities (or similar processes) provide a novel mechanism 
for heating the universe. Integrating over the energy flux per 
mean free path of all known TeV blazars yields a luminos- 
ity density, or equivalently a local heating rate, that dominates 
that of photoheating by more than an order of magnitude at 
the pr esent epoch, after a ccounting for incompleteness correc- 
tions dChang et al.H2012l hereafter Paper II). We have demon- 
strated that the local TeV blazar luminosity function is con- 
sistent with a scaled version of the quasar luminosity function 
dHopkins et al. 2007), thus the conservative assumption is that 
they evolve similarly, presumably due to the same underlying 
accretion physics. With this assumption, we showed in Paper 
II that for the redshifts at which blazar heating is likely to be 
important, z < 3.5, the heating rate will be relatively uniform 
throughout space. Between z ~ 3.5 and 6 it may experience 
order 50% fluctuations, and by z > 6, it will exhibit significant 
stochasticity with order unity deviations. 

This heating differs from other feedback prescriptions in an 
important way: since the number density of EBL photons and 
that of TeV blazars are nearly homogeneously distributed on 
cosmological scales, so is the resulting pair density. Hence, 
the implied heating rate is also homogeneous, i.e., the volu- 
metric blazar heating rate is uniform and independent of IGM 
densitjQ. Due to the large mean free path of TeV-photons, 
which is much larger than the turn-around radius of any virial- 
ized structures, the heating process is not expected to be domi- 
nated by contributions in highly biased regions at late times for 
z < 3.5. (At early times, when blazar heating exhibits consid- 
erable fluctuations, clustering bias could substantially modify 
the phenomenology of the heating mechanism.) The effect of 
a uniform heating rate is that the energy deposited per baryon 
is substantially larger in more tenuous regions of the universe. 
As a result, underdense regions experience a larger tempera- 
ture increase, producing an inverted temperature-density rela- 
tion in voids, asymptotically approaching T oc p~ l . Generally, 
we found in Paper II that without any fine tuning it is possi- 
ble to reproduce the inverted temperature-density relation at 

4 The term "blazar heating" exclusively denotes the dissipation of gamma- 
ray luminosity of HBL and IBL blazars between 100 GeV and 10 TeV with an 
equivalent isotropic luminosity of L > 10 42 erg s _I (see Paper II for details). 
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z = 2 . — 3 inferred by h igh-redshift Lya studies dBolton et al.l 
2008]; IViel et all 120091) . while simultaneous ly satisfying the 
temperature constraints at z = 2 (e.g., those bv lLidz et al.ll2010T) 
and leaving the local Lya forest unaffected. 

In a follow-up paper by lPuchwein et al.l (1201 ll) . we used hy- 
drodynamic simulations of cosmological structure formation 
to explicitly demonstrate that blazar heating provides not only 
an excellent description of the one- and two-point statistics of 
high-redshift Lya forest spectra but also the line width and col- 
umn density distribution. This detailed agreement includes re- 
producing the observed mean transmission, and is achieved us- 
ing the most recent estimate of the evolution of the photoion- 
izing background without any tuning. These successes are due 
specifically to the salient properties of blazar heating, in partic- 
ular its excess energy injection into the low density IGM and 
its continuous nature. 

In this work, we are interested in the impact of such an im- 
portant heating mechanism on cosmological structure forma- 
tion. While we propose a well-motivated, physical heating 
mechanism that uses a certain type of plasma instability, our 
conclusions concerning the thermal history of the IGM, the 
Lya forest, as well as structure formation remain robust in that 
they only rely on a mechanism that dissipates the energy of 
TeV blazars independently of the density and employs a red- 
shift evolution similar to that of the quasar luminosity func- 
tion. To study the impact on structure formation, we turn to 
the evolution of the entropy because in the absence of radia- 
tive cooling, entropy is conserved upon adiabatic compression 
and can only be increased through dissipation of gravitational 
energy in structure formation shocks or heating due to photo- 
ionization. Hence, entropy is a unique thermodynamic vari- 
able with which to learn about the impact of non-gravitational 
feedback processes. We identify two different classes of ob- 
jects where the time variable minimum entropy, or entropy 
floor, of the IGM induced by blazar heating might dramati- 
cally change our present picture of structure formation: the 
structure of galaxy groups and clusters, and the formation of 
dwarf galaxies. Before addressing the consequences of blazar 
heating for each, we introduce each separately, highlighting the 
relevant outstanding problems. 

1.1. The entropy problem in galaxy groups 

In the absence of non-gravitational energy injection, the X- 
ray luminosity of clusters (L x ) is expected to exhibit a self- 
similar scaling with intracluster medium (ICM) temperature 
(T) as L x ex T 2 from purely gravitational and shock dy namics 
jKaiserll986tlEvrard & Henry|1991l:lEvrard et aljl996l) . How- 
ever, this cannot extend from clusters down to groups without 
vastly overproducing the extragalactic soft X-ray background. 
Indeed, groups are observed to have a smaller X-ray luminos- 
ity compared to the self-similar exp ectation yielding a steeper 
scaling with temperature, L x oc T 3 (iMarkevitchl 1 1 9981) . This 
can be explained if the gas was initially preheated to an en- 
tropy floor of ~ 100 keVcm 2 whi ch reconci les the simulated 
background with observations (see Voit 2005, for a review). 

In principle, there are three physical processes that could 
produce such an entropy floor for groups, all of which may 
occur in practice, raising the question of which, if any, 
is dominant: (1) catastrophic cooling and collapse of the 
low-entropy gas at the centers of halos, allowing accretion- 
shock heated gas to adiabatically flow inward a nd replace the 
conden se d gas with its e levated entropy level dVoit & Brvanl 
12001 alibi : IVoit et all 120031) . (2) an early epoch of global en- 



tropy injection prior to the formation of gro ups and clus 
ters, typically referred to as "preheating" (e.g.. iKaiserl 1991 



Evrard & Henrvll99U|pWian et alJl999tlBalogh et al.ll 199 



9; 



Peril 1999HBorgani et alJl200lUBrvan & Voitll200 It ICroft et al 
20011: iBialek et al.l 120011: iBabul et al.l 120021: 1 Voit et alj I2003L 
2001 IBorgani & Viell I2009L IStaneketalJI2010l) . and (3) 
self-regulated AG N feedback a t the cores of groups 
and clusters (e.g., IChurazov et alj I200U ISiiacki & Springell 



2006t ISiiackiet al.l 120071 120081: IMcNamara & NulsenI 12007 



Puc hwein et al.l 120081: iBooth & Schavel 120091: iMcCarthy etafl 
2010tlDubois et al.ll2010tlTevssier et al.ll201 lb . 

Catastrophic cooling has been found to be unstable in large- 
scale numerical hydrodynamic simulations, typically result- 
ing in an untenably large fraction of baryons in stars (see 
Borg ani & Kravtsovl 120091 for a review). More importantly, 
runaway cooling generally precludes the existence of the ob- 
served hot groups, making it clear that it cannot be the sole 
explanation of the low X-ray luminosity of groups. However, 
radiative cooling is a physical process that will occur if it is not 
perfectly balanced by some heating process and the observa- 
tions of multi phase gas in cool core s (CCs) suggest that cooling 
is happening (Do nahue et al.H201 11). presumably through ther- 
mal instability dMcCourt et alj2012h . Hence, while some of the 
elevated entropies in groups have to be attributed to cooling, the 
absence of a cooling catastrophe implies a (non-gravitational) 
energy feedback process that will inevitably also raise the core 
entropy. 

In contrast, preheating has met with broad success in ex- 
plaining the observed L x - T correlation, ostensibly by sup- 
pressing the formation of dense cluster cores at low masses and 
therefore reducing L x for group-sized objects. Typically, the 
redshift range adopted for the injection of the entropy, z > 3, 
ensures that preheating is complete well before any galaxy 
group or cluster has turned around from the Hubble expansion 

-2/3 

and started to collapse. Hence the entropy, K e = kTn e ' (where 
n e is the electron density and k is the Boltzmann constant), is 
generally injected at the lowest possible gas density (and thus 
temperature) at that redshift. This minimizes the heating nec- 
essary to produce the observed impact upon the entropy, typi- 
cally 1 keV per particle. Numerical simulations incorporating 
a minimum entropy of the gas, or entropy floor, of K e ~ (100- 
200) keV cm 2 at redshifts around z — 3—4, are able to steepen 
the L—T relation from the purely gravitational collapse esti- 
mate (L t oc T 2 ), and find broad agreement with observations 
jKaisdll99ltlBalogh et d.lll999t iBiaiek et alj|2001b IVoit et al.1 
120021 120031: iStanek et al.ll2010l) . Nevertheless, there are seri- 
ous observational difficulties confronting the standard preheat- 
ing scenarios, apart from the theoretical problem of lacking a 
physical model for the heating process. First, after preheating, 
star formation i n and lower mass galaxies is significantly 
suppressed (e.g., IOh & Bensonll2003l: iBenson & Madaull2003b 
and hence inconsistent with observational data that finds, e.g., 
the peak of the star formation rate at redshifts z ~ 2. Sec- 
ond, the presence of a high entropy floor everywhere in the 
universe makes it impossible for groups to radiate away the ex- 
cess entropy and impossible to explain the existence of a sub- 
set of X-ray luminous, CC groups with steep entropy profiles , 
and low values of the central entropy (|McCarthv et al.l [2008; 
iFang & Haimanll200llCavagnolo et al.ll2009l) . 

Heating via feedback processes, e.g., from star formation 
(stellar winds and superno va) and AGNs (for a review, see 
IMcNamara & NulsenI 120071) . has already been shown to have 
a large impact on the formation and history of galaxy clus- 
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ters and groups. Self-regulated, inhomogeneous energy feed- 
back mechanisms are very successful in globally stabilizing 
the group and cluster atmospheres, and in particular, prevent- 
ing the cooling catastrophe. The resulting gas mass profiles, 
gas fractions, and t he L, - T correlation compare impressivel y 
to those observed (ISijacki et all 120081: IPuchwein et all 120081) . 
While there seems to be a globally convergent scenario emerg- 
ing with convincing energetics and duty cycles, the actual heat- 
ing process has yet to be identifie d and may involve in terest- 



nfilin et alJl201 lh . or plasma in- 



ing astro physics, e.g ., turbulence (Enfil in & Vogtl 120061) . cos- 
mic rays dGuo & Ohll2008b IT" 
stabilities dKunz et al.ll2011l) . 

As we will show in this paper, the entropy injection by 
blazars is in some sense an amalgam of both the preheating and 
feedback mechanisms. Unlike typical AGN feedback, the ef- 
fect of the blazar heating is not localized, operating at much 
larger distances typically and thus over considerably longer 
timescales. Unlike instantaneous preheating, the blazars pro- 
vide a time-dependent entropy injection rate, peaking near 
Z ~ 1 . As a consequence, we will show that the formation of 
L* galaxies is not suppressed and early-forming groups are lit- 
tle affected by blazar heating, having had time to cool and de- 
velop an X-ray luminous CC; avoiding the primary difficulties 
with the standard preheating scenario. On the other hand, the 
blazar-heated ICM ending up in late-forming groups will not 
have had sufficient time to radiate the additional entropy away 
before it can be gravitationally reprocessed in merging shocks 
which are ubiquitous in a hierarchically growing universe. We 
will argue that this can lead to elevated entropy core values re- 
sembling those of non-cool core (NCC) clusters. 

1.2. The dwarf problem in our Galaxy and nearby voids 

The A cold dark matter (ACDM) concordance cosmology 
predicts that Milky Way-sized halos should contain many more 
dwarf-sized halos than the observed number of dwarf galaxies, 
the so-called "substructure problem", or within the local con- 
text the "missin g satellites problem" (for a recent review, see 
iRravtsovl 120101) . Closely related to this problem is the "void 
phenomenon" which is the apparent discrepancy between the 
number of dwarfs in low-density regions in simu lations and the 
paucity of those in observations (Peebles 200 lh . In principle, 
both problems can be solved in three general ways: by sup- 
pressing the formation of dwarf halos, suppressing the forma- 
tion of galaxies within existing dwarf halos, and/or suppressing 
star formation within dwarf galaxies. In the case of the "void 
phenomenon", there is a forth class of models that solves the 
problem. Considering interacting dark matter and dark energy 
as mediated, e.g., by a Yukawa coupling, implies a fifth force 
that reduces late accretion onto halos and pushes matter out of 
voids, hence resolving the discrepancy of dwarf ab undances in 
voids dFarrar & Peeblesll2004b iNusser et all 120051) . In the fol- 
lowing, we review the three classes of models that are able to 
solve both problems simultaneously. 

Suppressing the formation of dwarf halos is difficult to 
accomplish within the context of ACDM, requiring modifi- 
cations to the standar d cosmological paradigm, such as in- 
teracting dark matter ( ISper gel & Steinhardt 2000), modifica- 
tions to the seed perturbation spe ctrum (IZentner & Bullock! 
20031). or warm dark matt er (WDM lDalcanton & Hoganll2001b 



Mac cio & Fontanojl2010l) . While we will discuss the last of 



these in Section 13.31 given the current success the ACDM 
mod el has had predicting the observed halo s tructures (see, 
e.g., iDalal & Kochanekl 120021: iMao et al.l 120041) , we will not 
comment upon these possibilities further. 



Suppressing the formation of dwarf galaxies, i.e., pre- 
venting the accretion of baryons by existing dwarf halos, 
may be accomplished in principle by a v ariety of mecha- 
nisms, including phot oionization heating (lEfstathioul 11992b 
Kauffmann et alTl! 993b lOuinn et all 1 19961 iThoul & Weinberg! 
1996b iKitavama & Ikeuchil l2000t iBullock et al.l |2000[ 12001b 
Chiu et al.H2001b ISomervilldl2002tlDiikstra et al.ll2004l) or ac- 



cretion shock heating ( Scannapieco et al.ll2001blKravtsov et al.l 



l2004t iSigward et al.ll2005l) .~ If baryons can collect, they may 
not be able to cool efficiently due to a l ack o f H I as a result 
of photoionizatio n dHaiman et al.1119971 120001) or intrinsically 
low metallicities dKravtsov et al.H2004b iKaufmann et al.ll2007b 



Tass is et al.ir2008t iRobertson & Kravtsovl T2008: Gn edin et al.l 
20091) . Nevertheless, recent numerical simulations which self- 



consistently include the photoionizing background due to star 
formation have found that while it does have a pronounced ef- 
fect, it cannot suppres s dwarf galaxy formation at the level im- 
plied by observations dHoeft et al.ll2006b lOkamoto et al.ll2008b 
iNickerson et al.ll201 ll) . More importantly, the metallicity dis- 
tribution of some dwarfs is i nconsistent with dwarfs generally 
being pre-reionization relics (iDolphin et al.ll2005blFenner et al.l 
12006b (Holtzman etafl 12006b lOrban et al.ll2008l) . and therefore 
the suppression of dwarf formation must have occurred more 
recently than the epoch of reionization. 

Once formed, the gas may b e removed from dwarf 
galaxies via photo -evaporation (Bar kana & Loebl [1999; 
IShapiro et al.l 120041) and feedb ac k fro m the first su- 
perno va e dMac Low & Ferraral fl999l: iDekel & Wool 



2003b iMashchenko et al.l 12008b Uube gas e al.l 12 008: 



Wadepuhl & Springelll201 lb INickerson et al. 
20121) . Tidal interactions of satellite dwarf halos with 



gas e al.l IZUUift 
l201lHUhligetaLl 



the Milky Way may result in a dramatic decrease in their 
mass, and to a lesser extent in circular vel ocity, after z ~ 2 
dKravtsov et al.ll2004c INickerson et al.ll201 ll) . However, these 
processes also strip material from larger halos, with the result 
that the smallest dwarf spheroidal galaxies presently within 
the Local Group may have had a mass at formation that was 
much larger than currently observed, and thus were capable 
of building up a sizeable stellar component in their seemingly 
shallow potential wells. 

As we will show in this work, the heating due to blazars 
provides an additional mechanism to suppress the formation 
of dwarfs. Unlike photoionization models, which typically in- 
voke the heating at reionization, blazar heating provides a well 
defined, time-dependent suppression mechanism, with the sup- 
pression rising dramatically after z ~ 2. In addition, due to 
its insensitivity to density, the heating from blazars suppresses 
structure formation most efficiently in the low-density regions 
that are responsible for late-forming dwarf halos. As a result, 
the impact from blazars is not degenerate with variations in the 
parameters of reionization and/or tidal interactions. 

1.3. Structure of this Paper 

This is the third in a series of three papers that discuss the po- 
tential cosmological impact of TeV emission from blazars. Pa- 
per I provides a plausible mechanism for the local dissipation 
of the TeV luminosity, effectively producing an additional heat- 
ing process within the IGM, and its implications upon high- 
energy gamma-ray observations. Paper II estimates the magni- 
tude of the new heating term, describes the associated modifi- 
cations to the thermal history of the IGM, and shows how this 
can explain some recent observations of the Lya forest. Pa- 
per III, this paper, considers the impact the new heating term 
will have upon the structure and statistics of galaxy clusters 
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and groups, and upon the ages and properties of dwarf galaxies 
throughout the universe, generally finding that blazar heating 
can help explain outstanding questions in both cases. 

In Section|2l we show that blazar heating necessarily implies 
the injection of a tremendous amount of entropy and employ 
the implication for the evolution of the warm and hot phases of 
the IGM. In particular, we show the broad implicatio ns of this 
heating on the X-ray population of clusters (Section [2. U , the 
entropy structure of galaxy groups/ clust ers (Section [2. 21 . the 
CC/NCC cluster bimodality (Section l2~3l ). and its im pact on the 
Sunyaev-Zel'dovich (SZ) power-spectrum (Section I2.41 i. We 
discuss the impact of this heating on the formation of structure 
in the universe in Section [3] showing from the viewpoint of 
linear theory that the formation o f dw arf galaxies that end up 
in Milky Way -sized halos (Section lXTT ) as well as those in voids 
(Section [3~2l are suppressed at late redshift (z < 2). Finally, we 
conclude in Section@] 

The c alculations presented below assume the WMAP7 cos- 
mology (lKomatsuetal.ll2011[) with h = 0.704, fi DM = 0.227, 
n B = 0.0456, ft A = 0.728, cr s = 0.81, n s = 0.967, and a mat- 
ter transfer function th at accounts for the baryonic features 
(lEisenstein&Hulll998l) . 

2. THE EVOLUTION OF ENTROPY AND THE HOT PHASE OF THE 
INTERGALACTIC MEDIUM 

The injection of large amounts of energy into the IGM by 
blazars is accompanied by a substantial increase in the IGM 
entropy. Here, we discuss the impact of this additional en- 
tropy upon the evolution of the gravitationally heated warm 
and hot components of the IGM. In doing so we neglect ra- 
diative cooling, an approximation that is well justified in the 
low-density regions of the IGM generally and one that we will 
justify for late-forming groups and clusters below. In this case 
the entropy of the IGM necessarily increases during gravita- 
tional collapse (e.g., due to structure-formation shocks, feed- 
back, and/or photoheating). As a consequence, the entropy in- 
jected by TeV blazars places an elevated floor upon the entropy 
of late-forming structures within the IGM. If radiative cooling 
is permitted (e.g., at the centers of early-forming clusters), the 
gas entropy may decrease^ reducing the impact of the blazar 
heating. In this work, we will encounter two different defini- 
tions of entropjQ which are proportional to each other, namely, 



K, 



kT 



and K: 



kT 



5/3 2/3 ' 

Ps m P Pg 



(i) 



Kg is a quantity conveniently used in the X-ray literature as it 
can be directly constructed from the observables temperature, 
T, and electron density, n e . The quantity K is the constant of 
proportionality in the equation of state P(p g ) for an adiabatic 
monatomic gas with mass density p g and pressure P. 

The evolution of the temperature and the entropy are shown 
in Figure Q] for a fluid element at overdensity 8 = p/ p- \ = 
(where p and p denote the average and local matter density) for 
the cases of pure photoheating and two realizations of blazar 

5 According to the second law of thermodynamics, the entropy of an isolated 
system which is not in equilibrium has to increase over time and can only 
remain constant during a reversible process. In the case of radiative cooling, the 
gas shares its entropy with the cosmic radiation field — a process that effectively 
removes entropy from the gaseous phase. However, the total entropy of the 
system including gas and radiation increases in this process. 

6 These definitions of entropy are related to the standard thermodynamic 
definition of entropy per particle by s = klnK 3 / 2 +so, where so is a constant 
that depends upon fundamental constants and a mixture of particle masses. 



heating. Our standard model is normalized to the observed 
number of TeV blazars, accounts for an incompleteness cor- 
rection for incomplete sky coverage at TeV energies and duty 
cycle, employs conservative assumptions of spectral correc- 
tions and contributing source classes for blazar heating, and 
has a heating rate of 7 x 10~ 8 eV cm~ 3 Gyr _1 at z = (Paper 
II). The heating rate in the optimistic model is a factor two 
larger, implying a rate of 1.4 x 10~ 7 eV cm~ 3 Gyr _1 at z = 
and matches the inverted temper ature-density relation found in 
high- r edshift Lya observatio ns dBolton et al.ll2008t IViel et all 
120091: iPuchwein et alJl20H . Since any mechanism that in- 
jects energy will produce a corresponding increase in K e , in 
the absence of cooling the entropy accumulates monotonically 
regardless of the mechanism responsible for heating the IGM. 
Nevertheless, some general statements can be made about the 
distinction between photoheating and that due to blazars. 

Due to the ionization balance maintained between recombi- 
nation and photoionization, photoheating produces a generic 
entropy injection profile following the e poch of reionization ; 
this is simply the "los s of memory" effect iHui & Gnedinl 19971 ; 
iHui & Haimanll2003h . In the absence of the blazar component, 
photoheating of hydrogen is the main source of heating for the 
universe at mean density at late times. The heating rate is given 
by the rate at which hydrogen recombines and is then reion- 
ized. In photoionization equilibrium, the temperature of the 
system will adjust such that the timescales of the net photoheat- 
ing and recombination cooling rates become equal implying a 
steady state. Since recombination does not depend on the UV 
photon density the rate of photoionization is also independent 
of the UV photon density. Hence the late time photoheating 
of the IGM is independent of the number of ionizing photons, 
and thus the number of ionizing sources (above the ionization 
threshold). As a consequence, the photoheating contribution to 
the IGM entropy at z = is limited to near K e pa 8keVcm 2 for 
6 = 0. 

Such an argument does not work for blazar heating. Here 
the efficiency is close to the maximum of 100% all the timefl 
provided there are sufficient numbers of EBL photons and each 
point in space is reached by a number of blazar beams of TeV 
photons (which is likely the case; see Section 3.2 of Paper II). 
Therefore, the heating rate, and thus temperature of the IGM, is 
nearly linearly dependent upon the VHEGR luminosity density 
of TeV blazars. Correspondingly, the entropy injection from 
blazars depends sensitively upon the history of the TeV blazar 
population, and thus the cumulative contribution of blazars to 
the present-day K e is somewhat uncertain. Nevertheless, given 
our conservative estimate of the blazar luminosity density in 
Paper I, fixing it to the quasar luminosity density, we find that 
blazar heating raises the entropy substantially, starting around 
the epoch of He II reionization (z ~ 3.5) and by z = 0, the in- 
clusion of blazar heating raises the entropy of the 6 = fluid 
element to K e ss (50- 100)keVcm 2 , approximately an order of 
magnitude larger than the case of photoheating alone. 

7 In general, the efficiency for heating depends on the saturated, nonlin- 
ear damping rate of the maximally growing mode of the dissipating plasma 
instability, Tm, and the inverse Compton cooling rate, Tic, and scales as 
FtviArM + Tic). We assume that the nonlinear damping rate is equal (or of 
order) the linear growth rate (see Section 3.5 of Paper I). For the linear growth 
rate of the oblique instability, this efficiency is close to unity for luminous 
blazars (Paper I). Numerical simulations for a mildly relativistic pair beam pen- 
etrating into a hot, dense background plasma suggest that a significant fraction 
(~ 20%) of the beam energy is heating the background plasma through the 
oblique instability before the two-stream plasma instability takes over, poten - 
tially dissipating another fraction of the beam kinetic energy (Bret et al. 2010). 
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Figure 1. Thermal history (left) and entropy history (right) of a patch of mean density (<5 = 0) of the IGM. The solid curves are for pure photoheating with sudden 
reionization histories for H and He II at i lc ion = 10 and zhcii = 3.5. The effect of the "loss of memory" in the thermal history of the IGM is apparent in the entropy 
as the pure photoheating curve asymptotes to K e ss 8 keV cm 2 at z, = 0. The dashed (dash-dotted) lines show the evolutio n for the standard (op timistic) blazar heating 
model that employs the blazar luminosity density, i.e., using the redshift evolution of the quasar luminosity density from Hopkins et al. (2003) an d are normalized to 
the local heating rate. At z = 0, those differ by a factor of two and the standard one has been estimated from the observed number of blazars, emplo ying a conservative 
estimate for the incompleteness correc tion, while the optimistic one naturally matches the inverted temperature-density relation in the Lya data I Bolton et al. 2008; 
IVieletal.120 09: Puchwein et al. 20 lj). 



• z=0.5 




0.1 1 10 

1+5 



Figure 2. Entropy-density relation (Ke versus 1+5) for z = 0.5 (red diamonds), 
1 (green triangles), 2 (blue squares), 3 (black circles), and 4 (cyan inverted 
triangles) for a pure photoheating model. 

As discussed in Paper II, the magnitude of the temperature 
enhancement due to blazar heating is density dependent: while 
the photoheating rate depends linearly on density, blazar heat- 
ing is independent of density and dominates photoheating for 
regions with 5 < 10. Thus, blazar heating is expected to have 
the largest effects in voids. In Figures [2] and Q] we show the 
entropy-density relation (K e vs. 1 +5) for a variety of redshifts, 
ranging from z = 0.5 to z = 4. For comparison, Figure [2] shows 
the case of photoheating only: the entropy-density relation re- 
flects the effect of photoionization which deposits a uniform 
energy per baryon during H and He II reionization implying 
K oc (1 +<5)~ 2 / 3 . This effect is quickly erased due to the "loss 
of memory" effect as well as adiabatic cooling due to the Hub- 
ble expansion approaching an asymptotically constant entropy- 
density relation at late times. 



Figure [3] shows the effect of TeV blazar heating on the 
entropy-density relation for our standard (top) and optimistic 
(bottom) model. As expected, blazars have the most dramatic 
effects upon K e in low-density regions (1+<5~0.1) due to the 
spatially uniform and IGM density independence of their vol- 
umetric heating rate. This steepens the entropy-density rela- 
tion, causing it to approach an asymptotics of K oc (1 + <5)~ 5//3 
in voids. In fact, blazar-induced entropies can exceed those 
due to photoionization by a factor of 125-250 in these low- 
density regions (1 + 5 ~ 0.1, see right-hand panels of Figure 
0, reaching K e ~ (1250-2500)keVcm 2 by z = 0.5. For high 
density regions (1 + 5 ~ 10) the effect is still pronounced, in- 
creasing the IGM entropy by up to an order of magnitude, i.e., 
K e « (10-20)keVcm 2 at z = 0.5. 

Apart from increasing the mean entropy, blazar heating also 
increases the scatter in entropy at fixed density which is es- 
pecially noticeable for z = 0.5. Fractionally, this scatter ap- 
pears to be roughly twice as large when blazar heating is in- 
cluded, though in absolute terms the scatter is much larger for 
the blazar cases. This is because even in linear theory, the en- 
tropy (and temperature) of any patch of overdensity 5 depends 
on its "collapse" history. Namely, whether it collapses like a 
Zel'dovich pancake, i.e., in one dimension first which implies 
a smaller entropy and temperature, or more spherically which 
produces a greater entropy and temperature, i.e., more akin to 
accretion onto a local overdensity. 

The fractional scatter in Xbiazar/^photo is much larger than 
that found in either entropy individually. The reason is that the 
scatter induced by blazar heating is anti-correlated with that in- 
duced by photoheating due to their very different dependencies 
upon IGM density. Patches that begin at low densities when 
blazar heating is ignored are biased toward lower entropies 
as a consequence of an extended period of slow recombina- 
tion. When blazar heating is included, these same patches are 
more efficiently heated, and therefore biased toward higher en- 
tropies. 

The redshift evolution of blazar heating introduces additional 
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Figure 3. Entropy-density relation (K e versus 1 + 6, left) for the standard (top) and optimistic (bottom) blazar heating model. Blazar heating increases the entropy 
over the pure photoheating case by a factor ^blazar/^photo which is shown versus 1+5 (right) again for the standard (top) and optimistic (bottom) blazar model. 
Shown are different redshifts z = 0.5 (red diamonds), 1 (green triangles), 2 (blue squares), 3 (black circles), and 4 (cyan inverted triangles). 



stochasticity: patches at a given density sample a distribution 
of turnaround times and hence preheated entropy values. Af- 
ter turnaround, the gas is adiabatically compressed and moves 
on an adiabat (a line of constant K e ) to higher density with 
an entropy value that depends on the "collapse" time, and the 
steeper entropy-density relation generates the larger scatter at 
any 5 > 1. Hence by preheating the universe, the first shells of 
gas that are able to collapse onto an object do not experience 
shock heating as they are only adiabatically compressed which 
results in an entropy floor of the object after formation. The 
later collapsing shells experience weaker shocks with a smaller 
Mach number due to their already elevated entropy level. 

Associated with the large increase in the entropy of the IGM 
are a number of observational effects. In the following we dis- 
cuss the impact blazar heating has upon the correlation between 
cluster X-ray luminosity and temperature, the entropy profiles 
and the CC/NCC bimodality of clusters and groups, and the SZ 
power spectrum. In each we note the unique effects imposed 
by the relatively recent nature of the entropy injection. 



Blazar heating raises the entropy floor for z < 2 with its 
peak contribution around z ~ 1 and acts on scales much larger 
than the turnaround region of clusters. The formation epoch of 
groups and clusters is roughly matched to the epoch of heating. 
As a consequence, early-forming, i.e., old groups may be lit- 
tle affected by blazar heating, had time to cool and develop an 
X-ray luminous core, potentially representing the class of CC 
groups/clusters. On the other hand, the blazar-heated IGM col- 
lapsing into late-forming, young groups will not have had suffi- 
cient time to cool. These groups will "remember" the elevated 
entropy floor as an extended core that substantially changes the 
initial conditions for their subsequent hierarchical evolution. 
However, because the heating occurs at late times, and thus af- 
ter the first groups and clusters have begun to form, a larger 
energy input is required in comparison to the traditionally em- 
ployed early preheating models. 

For the observability of such a non-gravitational entropy 
floor and its implication on the thermal history of the ICM, it is 
important to consider its cooling timescale for typical entropy 
values of the IGM around z — 0.5 and parameters typical of the 
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intra-group medium: 



^cool 



3nkT ^ ( K e 

2n e ii H A(T,Z) ' yr X V75keVcm 2 



3/2 



IcT 
IkeV 



-1/2 



A(AT,Z) 



A(lkeV,0.3Z Q ) 



(2) 



Here n = 2.2«h, A(T,Z) is the cooling function at a given T 
and metallicity, Z, (whi ch is a relatively flat function of T near 
AT = 1 keV and 0.3Z Q , ISutherland & Dopita|[l993l) . A look- 
back time of 4.5 Gyr corresponds to a redshift of z = 0.425. 

Instantaneous preheating by blazars (see Figure [TJ is in- 
capable of directly producing the highest entropies (up to 
600keVcm 2 ) that are presently observed in a few clusters 
(ICavagnolo et al.l [20091) . Comparing the cooling time to the 
dynamical time of a cluster/group of approximately ~ 1 Gyr or 
the typical timescales of significant turbulent pressure support 
of ~ 2 Gyr after a cluster merger w hich provides continuous 
heating throughout dPaul et al.ll201 ll) . we expect our heating to 
indirectly impact the entropy distribution through gravitational 
reprocessing of the blazar preheated core entropy, leading to a 
gravitational amplification of entropy that we now explain. 

A larger central entropy, or equivalently, a lower central den- 
sity of a merging cluster facilitates shock heating which im- 
plies an increase of the core entropy of the final object com- 
pared to that of a dense cooled core. CCs have a large inertia 
causing them to be rather resilient against shock heating, typi- 
cally surviving the merger with only a marginally elevated en- 
tropy level which is then subject to fast radiative cooling. This 
effect gives rise to the well-known overcooling problem in cos- 
mological simulations of galaxy clusters i f cooling is not coun- 
teract ed by any feedback process (e.g., iBorgani & Kravtsovl 
12001 . 

In order to quantify the effect of gravitational reprocessing of 
a preheated entropy core, let us define a net entropy amplifica- 
tion factor through shock heating. This is easiest done by com- 
paring two types of simulations of galaxy cluster formation: 
one with gravitational physics only and one where gravitational 
physics is supplemented by a preheating epoch of the IGM at 
mean density with a uniform entropy floor of Kf\ oor which pre- 
cedes the turnaround and formation of a galaxy cluster. In the 
first case, the core entropy structure is set by gravitational for- 
mation shocks leading to a central entropy of Xg rav .o- In the 
preheating case, the core entropy is raised to an elevated level 
after nonlinear evolution during virialization, Ap, e .o- Hence, the 
net entropy amplification factor is then defined as the ratio of 
the core entropy in the preheating case, K pie x), to that obtained 
by adding the IGM entropy floor value to the entropy in gravi- 
tational heating case only, Kf[ oor + K gmv ,o- 

Non-radiative cosmological simulations of galaxy clusters 
and groups demonstrate that the net entropy amplification fac- 
tor can reach values ranging from 3 to 5 for clusters and groups, 
respectively, in the high-entropy case with a preheating value 
of -^floor = lOOkeVcm 2 imposed at z = 3 (see Figures 1 and 2 
in IBorgani et al.ll2005l) . For the low-entropy case with Kfi oor = 
25keVcm 2 , this amplification factor is only mildly reduced to 
2-5 for clusters and groups, respectively. This entropy ampli- 
fication factor seems to be reduced for radiative cluster simu- 
lations that include cooling and star formation, possibly due to 

8 For a cluster with a given mass, a larger core entropy implies a lower gas 
density as the (constant) core temperature kT oc Kn 2 / 3 reflects the virial value 
imparted by the cluster's potential depth. 
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Figure 4. Median mass accretion histories (top) and mass accretion rates de- 
fined by the logarithmic slope S = dlogM2QQ c /d\oga(z ) (bottom) for five char- 
acteristic halo masses Mjqoc according to the model by Zhao et al. 12009). Ar- 
rows label the half-mass redshifts of the systems. At late times, the mass accre- 
tion rates on group scales start to flatten out which is a result of the change from 
the fast to the slow accretion mode. In contrast, clusters are still forming today 
as the half-mass redshift zo 5 are still in the fast accretion regime characterized 
by S > 2. 

the assumed early epoch of instantaneous entropy injection at 
z = 3 that facilitate s cooling of the entr opy floor thereafter (see 
Figures 3 and 4 in Bo rgani et al.ll2~005l) . Groups, however, still 
remain severely affected by an early epoch of preheating. Thus 
in combination with gravitational amplification, blazar heating 
could have an important effect on the subsequent thermody- 
namic evolution of late-forming groups and possibly the clus- 
ters they evolve into. This, however, is subject to nonlinear 
structure formation and needs to be carefully studied with cos- 
mological hydrodynamic simulations containing a large sample 
of well-resolved galaxy clusters which sample the full distribu- 
tion of formation redshifts. 

The impact of the entropy injection due to blazar heating 
upon clusters is somewhat less clear. Clusters form in highly 
biased regions through the mergers of smaller, virialized sys- 
tems such as groups. In most cases, these need to have already 
formed by z = 1 which can be inferred from the mass accre- 
tion history in Figure [4] This shows that on average, the most 
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massive group progenitor of a cluster with virial mass3M2oo e = 
10 14 Mq crossed a mass threshold of 3 x 10 13 M before z = 1 . 
At this time the entropy floor due to blazar heating was still 
smaller with typical values K(5 = Q,z = 1) ~ (25-50)keVcm 2 
implying a smaller entropy core in groups at higher redshifts 
Z > 1. Hence, we naively may not expect blazar heating to have 
a large impact upon the X-ray luminosity and entropy profiles 
of clusters. However, this is not the case for two reasons. 

First, the mass accretion rate, defined by the logarithmic 
slope S = dlogM200c/d\oga(z), is larger for more massive sys- 
tems and at higher redshifts (where we have introduced the cos- 
mic scale factor, a = 1/(1+ z); see bottom panel of Figure |4|. 
The mass accretion history of a ha lo observed at z = grows 
on average as M(z) = M exp(-£z) (IWechsler et al.ll2002l) . The 
single free parameter in the model, £, can be related to a charac- 
teristic formation redshift 1 +z c = £/S- Hence we find ourselves 
in the fast accretion regime if the mass accretion rate is larger 
than a characteristic value usually taken to be S = 2. In partic- 
ular, the relative mass accretion rates increase from z = to 2 
by a factor 3 for clusters (M200C = 10 15 M Q ) and 10 for groups 
(M 2 ooc = 10 13 M o ) (see also iGottlober et al.ll200Tb . A larger 
relative mass accretion rate implies that the gas in the core is 
heated by accretion processes at a faster rate. This means that a 
lower preheated core entropy value early-on can in principle be 
gravitationally processed at a substantially higher rate and, as a 
result of this, move onto a higher adiabat with a longer cooling 
time. 

Second, blazars like AGNs follow the clustering bias of the 
matter density field and hence turn on first in highly biased re- 
gions, i.e., regions that evolve into clusters and super-clusters. 
At late times (z < 3.5) blazar heating is expected to be nearly 
spatially uniform, since more than a single blazar contribute 
significantly to the local heating rate of any given patch of the 
universe (Paper II). At much lower redshifts, the number of 
contributing sources grows dramatically, approaching 10 3 at 
z = 0, and thus today blazar heating is nearly homogeneous. 
Prior to z — 3.5, however, blazar heating may exhibit ~ 50% 
fluctuations locally, due to the paucity of blazars in the early 
universe. By z — 6 as much as 75% of the local heating can be 
due to a single object, implying large Poisson fluctuations in 
the heating rate. This results in a clustering bias at early times. 

The rare first blazars are expected to appear in highly biased 
regions, corresponding to those that later evolve into groups 
and clusters, and therefore themselves be clustered. The likeli- 
hood of a patch of the universe being heated by such a blazar 
is then enhanced in the biased regions, where both the heating 
rate is larger (since the VHEGR density is larger) and the prob- 
ability of being covered by a blazar in the cluster is largerFI 
It is this early preheating in groups that may facilitate their 
evolution into NCC clusters (via the gravitational reprocess- 
ing of the high-entropy cores). Therefore, while we will focus 

9 In this Section, we define the virial mass M200C as the mass of a sphere 
enclosing a mean density that is 200 times the critical density of the universe. 

10 In the limit of a few, highly clustered blazars, the VHEGR intensity de- 
clines exponentially from the clustering site for a spatially constant distribu- 
tion of EBL photons, giving rise to an exponentially decreasing heating rate 
(on the scale of the VHEGR mean free path which is larger than the clustering 
length scale). However, if the distribution of EBL photons is also clustered 
and correlates spatially with the highly biased regions (which is expected in a 
hierarchically growing universe), the heating rate of those biased regions will 
be additionally enhanced in comparison to the void regions that have a low 
probability of being covered by a blazar. We stress that the effect of spatially 
and temporarily inhomogeneous blazar heating due to the clustering bias is ex- 
pected to be absent at late times (z < 3.5) due to the frequent occurrence of 
blazars. 



upon group/small-cluster mass scale in our study of the entropy 
structure immediately after group/cluster formation, gravita- 
tional reprocessing combined with the bias in the early blazars 
suggest that these effects may be important for more massive 
clusters as well, propagating the effect of blazar heating up the 
mass hierarchy. 

The formation time of groups and clusters determines the 
blazar contribution to their entropy, inducing an intrinsic scat- 
ter in the L x -T relation associated with the distribution of col- 
lapse histories. As a result, we would expect to find system- 
atically higher entropies in younger groups and possibly clus- 
ters provided gravitational reprocessing was effective. There 
is some evidence that this is the case. Optically selected, 
and therefore young, group and cluster samples have on av- 
erage lower X-ray surface brightness and sm aller gas mass 
fracti o ns compared to X - ray-selected sam ples dMahdavi et alj 
l2000t iHicks et alj 12001 iDai et alj |20T(I . Optically bright, 
small to moderately massive clusters (kT > 4keV) at redshifts 
z = 0.6-1.1 are under-luminous in X-rays for a given T, which 
implies a reduced gas density and by exte nsion an enhance d 
core entropy by roughly a factor of two ( IHicks et "ail 120081) . 
Similarly, an X-ray stacking analysis of ROSAT data based 
on optically selected groups at low redshift (0.5 keV < kT < 
2keV) from the Two Micron All Sky Survey catalog finds 
systematically lower gas mass fractions than expected, / gas , 
(within an over-density of 500 times the cri tical density of the 
universe) and flatter temperature profiles dDai et all ^OlO fH, 
again implying larger entropies. Taken at face value, this 
is in conflict wi th careful X-ray studies using Chandra data 
dSun et al.ll2009l) if both samples are believed to represent the 
same underlying distribution. Alternatively, this is perfectly 
consistent if the entropy of clusters and groups varies with for- 
mation time, with the peak entropy injection occurring near 
z ~ 1 . Both are well matched to the properties of AGN feed- 
back generally, and blazar heating especially. 

2.2. Thermodynamic structure of galaxy groups and clusters 

The large injections of entropy necessarily influence the for- 
mation of cosmic structure. While the incorporation of an in- 
stantaneous entropy floor in numerical simulations does repro- 
duce the L x - T scaling relations well, it is ad hoc and fails 
to produce the observed entropy profiles in a subset of X-ray 
luminous, cool core groups with steep entrop y pr ofiles and 
low values of the central entropy (see Section 11. It . Thus, it 
is clear that the effects of heating the ICM are more compli- 
cated than the introduction of a constant, glob al entropy floor 
dMcCarfhv et alj 120081: IFang & Ha iman 2008b. Here we ex- 
plore whether blazar heating may be capable of producing the 
observed entropy profiles as well as some of the features of the 
distribution of large-scale structures observed. 

While large-scale numerical simulations are required to 
study the entropy evolution of groups and clusters in detail, 
we can use the conservative property of entropy to estimate 
the effect of blazar heating on the thermodynamics of galaxy 
groups and (to a lesser extent) clusters. Smooth accretion, in 
which cold gas enters the cluster through a spherically sym- 
metric accretion shock, results in self-similar entropy profiles 

1 1 We caution the reader that such a stacking analysis could have potential 
biases: the stacked X-ray spectrum in a given richness bin might be domi- 
nated by a few hot systems and fitting an average spectrum could then bias the 
temperatures and hydrostatic masses high and hence the gas fractions low. Ed- 
dington bias from more numerous smaller systems can additionally lower the 
average gas fraction at a given optical richness. Before drawing far reaching 
conclusions, careful mock analyses are needed to confirm these results. 
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( Voit 2005). Characteristic values for the physical parameters 
are given by their average values within the virial radius, 7?2oo 
(that we define as the radius of a sphere enclosing a mean den- 
sity that is 200 times the critical density of the universe): 

n e 200 = 2Q0x e X n f b — = 1.7 x 10 _4 £ 2 (z)cm- 3 
ftp 

kT2a) = ^^ = ^ [l0GH(z)M2QQc ^ 



2/?200 



= 11^1] keV 



2/3 



200 = 



6 x 10 13 M Q 

5/3 / 9 

kT 2 oo _ pm p < I 4ir G z M 2 ooc 



(3) 



2/3 



„2/3 oo 2 \15 (I +X H )f b H(z) 

^(M2Wc£Hz}\ 2l \ .. 2 

= 326 ke V cm , 

V6xl0 13 M Q y 

where ft, = fib/ £l m is the universal baryon fraction, the electron 
fraction is defined as the ratio of electron and hydrogen num- 
ber densities, x e = n e /rin = (Xn + 1 )/(2 Xu) = 1 . 1 58, in which we 
assumed a fully ionized fluid with a primordial hydrogen mass 
fraction X^ = 0.76, and p a - = p a {z) = 3H 2 (z)/{8irG) is the crit- 
ical mass density, in which the Hubble function H(z) is given 
by 



H 2 



: e\z) = (i +zfn m +(i +z) 2 (i - n m - n A )+ n A . (4) 



Of particular importance here is that the entropy scales as 
K e oc r 11 (Voit 2005). The predictio ns of this simple mode l 
agree well with numerical simulations dTozzi & Normanl2 001) 
and the entropy profiles in the outer regions of clusters inferred 
from recent X-ray obse r vations by Chand ra and XMM-Newton 
(ICavagnolo et alJl2009t iPratt et al.ll2010t) . Within the smooth 
accretion model the stratified profile is built from the inside 
out, with later accreted shells containing larger entropy due 
to gravitational heating at accretion shocks. Thus, if radia- 
tive cooling can be neglected, the entropy distribution of the 
gas in the cluster core reflects the entropy of the IGM imme- 
diately prior to the initial collapse of the group/cluster. This 
assumption is reasonable as long as the radiative cooling time 
(Equation (fJI) is longer than the time interval between succes- 
sive mergers, or eq uiva lently the mass accretion timescale (see 
discussion in Sect. 12. U . 

In practice, the smooth-accretion picture is over-simplified. 
Most of the ultimately accreted gas is not smoothly distributed, 
but rather contained within virialized substructures, and there- 
fore has already been shock heated, fundamentally altering the 
way in which the entropy profile is generated. This results 
in a more complex morphology of the dissipating structure- 
formation shocks, which exhibit a rich network of shock fronts 
at which the gravi t ational energy of the gas is dissipate d 
(iMiniati et alJ l2000t iRyu et alJ 120031: IPfrommer et"aT] 120061) . 
Nevertheless, despite this apparent chaos, non-radiative galaxy 
cluster simulations that have sufficiently high resolution find 
approximately self-similar entropy, density, and temperature 
structures outside of the core region ; independent of t he numer- 
ical method used dFrenk et al.ll 1 9991 : 1 Voit et al.ll2.Q05h . yielding 
a universal entropy profile of 



A;(r,z)=1.45A e ,2oo(z)(r//f 2 oo) 



1.2 



(5) 



ranges itself such that it finds itself at the bottom of the gravita- 
tional potential and eventually mixes with the surrounding halo 
gas as can be seen by rewriting the hydrostatic equation, 



dP 
dr 



GM«r) { P^ 15 



(6) 



The reason for this is that the low entropy/high density gas ar- 



Magnetic draping of cluster fields during gravitational set- 
fling provides a t hermal insulation of these low-entropy parcels 
(lLvutiko\il2006t iDursi & Pfrommer!l20pl IPfrommer & Dursil 
120101) . and hence this settling occurs adiabatically, producing a 
stratified core in which the hot, high-entropy halo gas remains 
at large radii. 

If the universe is preheated, the central entropy of a newly 
formed group or cluster is replaced with a flattened core. Ini- 
tially, the entropy of the core is set by that of the IGM at the 
time the second gas shell is accreted, forming an accretion 
shock and adiabatically compressing the gas in the first. Em- 
ploying only heating by blazars and structure formation shocks, 
we show the resulting entropy profiles for groups of a variety of 
virial masses and formation redshifts in Figure|5] exploring the 
distribution around the median object mass. These are based 
upon Equations (0 and © as well as the floor values implied 
by Figure Q] Interestingly, the blazar-induced entropy floor 
at z = 0.5 is comparable to the gravitationally established en- 
tropies in groups of 3 x 10 13 M Q at radii (0.2-0.4) 7?2oo for our 
standard and optimistic blazar models, respectively. The im- 
plied core sizes are similar to observationally accessible radii 
of ^7,-joo ~ 0.3 R 200 of optically selected clusters (see Figure 1 
of lPratt et alj |2010l). In our picture, these correspond to young 
clusters with recent formation times. 

Thus far we have implicitly assumed the instantaneous for- 
mation approximation, i.e., that we may identify a particular 
redshift at which to calculate K e ,200 and the core entropy level 
due to blazar heating. This is naturally identified with the red- 
shift at which the group/cluster has assembled half of its mass, 
after which strong structure-formation shocks develop. Figure 
2] shows the accretion histories for groups/clusters in the mass 
range of relevance here, noting explicitly the half-mass red- 
shifts for a variety of accretion histories. Typically, these occur 
after z — 1, by which time blazar heating has had an oppor- 
tunity to inject significant amounts of entropy into the IGM, 
implying that the TeV blazars can have a substantial impact 
upon the structure of groups/clusters in practice [3 However, 
we note that median mass accretion histories are too simpli- 
fied to assess the impact of blazar heating upon group/cluster 
entropy profiles in detail since they do not address the radial re- 
distribution of accreted gas. Furthermore, we do not attempt to 
address the generation of entropy by gravitational reprocessing 
during mergers of later accreted material in Figure|5] After this 
explicit demonstration of the impact of a redshift-dependent 
entropy floor on cluster entropy profiles, we turn to the impli- 
cations of these for the cluster population as a whole. 

2.3. Implications for the bimodality of core entropy values 

The thermodynamic properties of clusters in their centers 
show a clear bimodality, which is traditionally separated into 

12 The half-mass redshifts depend on the adopted definition for the halo 
mass (and of course cosmology which we fix here for simplicity). While 
we choose to use Mmoc in this section, we note that the half-mass red- 
shift are somewhat smaller when adopting the definition for the virial mass 
with an overdensity that varies with redshift ( Bryan & Norman 1998). In 
this case, we obtain zo.5 = {0.85,0.74,0.62,0.48,0.36} for our halos, M vir = 
{1.3,3.8,13,38,130} X 1O 13 M . 
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Figure 5. Radial profile of the entropy, K e = kT/nl' 3 , in galaxy groups/clusters for varying formation redshifts (left) and cluster masses (right) when our standard 
(top) and optimistic (bottom) blazar heating models are employed. While the outer profile results from dissipating gravitational energy in accretion shocks (dotted), 
the entropy core immediately after formation is set by the redshift-dependent blazar heating. Note that only late-forming groups (z < 1) are directly affected by 
this homogeneous but redshift-dependent preheating mechanism. If the ratio of cooling time, t mo \ to the time period to the successive merger fmeiger, is smaller 
than unity, the group will radiate away the energy associated with the elevated entropy core and evolve into a CC. Alternatively, if f coo ] > fmerger, merger shocks can 
gravitationally reprocess the entropy cores and amplify them (not shown here). Potentially those groups may evolve into NCC systems. 



two classes: CC and NCC clusters. The former are defined 
to have temperature profiles which decline significantly to- 
ward the center whereas the central temperature distribution 
of the latter remains constant and often correlates with merger 
events. As a result, the distribution of core entropy also appears 
to be bimodal in clusters, with the CC population peaking at 
K e fi ~ 15keVcm 2 and NCCs at K e , Q ~ 150keVcm 2 separated 
by a gap between K e o ~ (30-50)keVcm 2 . Roughly half of the 
entire population of galaxy clusters in the Chandra archival 
sample show cooling times that are lo nger than 2 Gyr and hav e 
high-entropy cores typical of NCCs ( iCavagnolo et al.l [20091) . 
This CC/NCC bimodality appears to be real and not due to 
archival bias as a complementary appro ach has shown with a 
statistical sample (S anderso n et al. 2009|). 

Before we discuss how blazars may impact on the relative 
abundance of cluster populations, we will shortly review how 
the currently favored hypothesis of AGN feedback compares 
to data. The possibility that AGN feedback can raise the core 
entropy K e ,o to values representa tive of NCCs on the buoyancy 
timescale (e.g., Guo & Oh 2009) is not supported by observa- 



tions. This is explicitly shown in Figure|6]which correlates K e> o 
with the volume work done by the expanding bubbles, E c . dw , and 
with the cavity power, P cav , in systems that show X-ray cavities 
inflated by AGN bubbles. The total energy required to inflate a 
cavity is equal to its enthalpy, given by 

E c , v = ^-PV tot = 4PV tot (7) 
7-1 

assuming a relativistic equation of state 7 = 4/3 within the bub- 
bles. The cavity power is estimated using P cav = £ C av Abuoyancy 
dBirzan et al.l2004URafferty et al.l2006l) Fl We compare the en- 
ergy used to inflate the cavities to the gas binding energy of the 



The buoyancy timescale could be either an overestimate of the true age 
since the cavity is expected to move outward supersonically during the early, 
momentum-dominated phase of the jet. It could also be an underestimate of 
the true age since magnetic draping provides an addit ional drag force slowing 
down the rise of the bubble fDurs i & Pfrommer 2QQH). 
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Figure 6. Correlating the observed core entropy values of clusters, K e .o, (taken from the gradient extrapolation method of Cavaenolo et al. 2009) with the AGN 
cavity energy £cav, 

as inferred from the volume work done by the expanding bubbles (left), and with the cavity power P cav = -Ecav / fbuoyancy (right, from lRaffertv et al.1 
120061) . Color coding reflects average X-ray temperatures, the lower limit of each color bin is labeled in the upper right of the left panel. Arrows denote the gas binding 
energy contained within a spherical region of radius R2500 — P200 1 3 ■ These very powerful AGN outbursts with energies up to 1 62 erg and powers of 1 46 erg s -1 are 
in some cases energetically capable of unbinding the gas in the core regions. However, the mechanical energy of these expanding cavities only heat the cluster core 
enough to prevent a cooling catastrophe. On the buoyancy timescale, no AGN outburst transforms a CC to a NCC cluster. This is apparent from the low core entropy 
values (with a median K e = 15keVcm 2 ) of typical CC clusters. 



core region within a spherical region of radius P2500 — ^20o/3, 

GM\ 5m = /g^ M 5/3 )[10G//(z)]2 /3 

(8) 



££,2500 = /gas,2500 



27? 



2500 



1 x 10 60 erg 



kT x \ 
IkeV / 



2 

3.23 



at z = 0. Here we use the phenomenological scalings obtained 
by X-ray observations of h 70 M 2 5oo = M s (fcr x /5key) L64± 006 , 
with M s = (2.5 ±0.1) x 10 14 /i^ M Q dVikhlinin et alj2006h and 
f ras .2s oo = (0.0347 ± 0.0016) (/tr x /lkeV) a509±0 034 (ISun et alJ 

P cav or £ cav both measure the energy and power that is in 
principle available for heating the ICM, if it can be efficiently 
tapped by some processQ However, as shown in Figure |6l 
even the most energetic and powerful AGN outbursts, with 
£ ca v ~ 10 62 erg andP cav - 10 46 ergs _1 , e.g., MS0735+7421 and 
Zwicky 2701, which are energetically capable of unbinding the 
gas in the core regions, are unable to disrupt the CC and trans- 
form the cluster into an NCC state on a buoyancy timescale. 
This is apparent from the low core entropy values (with a me- 
dian K e u = 15keVcm 2 ) of typical CC clusters. This is quanti- 
fied with a linear Pearson correlation coefficient of 0.71 corre- 
lating log K e x) with logP cav or \ogE cdv , respectively. This proves 
that the PdV work done by these expanding cavities is trans- 
ferred inefficiently to the surrounding medium, i.e., the ICM 

14 We note that of this 4PVtot, only PVtot is available in form of mechanical 
energy while the internal energy U = 3PV t ot is presumably still stored within 
the bubbles. What fraction of this internal energy is eventually thermalized, 
and thus can potentially contribute to unbinding the cluster gas and/or rais- 
ing the core entropy, is not a priori clear. If the energy is stored in cosmic 
rays or magnetic fields, it may be transferred to t he thermal pool via cosmic 
ray Alfven- wave heating I Kulsrud & Pearce 1 969) or magnetic reconnection, 
respectively. Nevertheless, the uncertainty induced by this is small in com- 
parison to the many orders of magnitude increase in cavity energy required to 
explain the non-cool core clusters. 



entropy produced by AGN inflated bubbles is far less that the 
virial value of K e ~ 540keVcm 2 x (M200C/M14) 2 / 3 , but enough 
to arrest overcooling. 

On the other hand, there is a strong anticorrelation between 
the radio power of the brightest cluster galaxy and K e o for 
nearby clusters (z < 0.2), implying that bright radio emis- 
sion is_greferentiafly "oji" for K e ,o < 40keVcm 2 (see Figure 
2 in iCavagnolo et all 120081) . While AGN feedback seems to 
be unable to transform a CC to an NCC cluster (on a buoy- 
ancy timescale), it appears to be critical in stabilizing the ther- 
mal atmospheres from entering a cooling catastrophe and col- 
lapsing. In principle the impact of AGN-induced turbulence 
on heat transport (conductively and advectively) could result 
in a CC t o NCC metamorphosis on a much longer (> Gyr) 
timescale dParrish et al.ll2010t iRuszkowski & Ohll2010f) . This 
is because the temperature difference between the maximum 
of the temp erature profiles to the cold center is at most a fac- 
tor of three ( IVikhlinin et alj |2006), implying that heat transport 
could initially increase the central entropy by a similar factor 
as Kg Q oc kT. In the absence of radiative cooling, the asso- 
ciated pressure enhancement would adiabatically expand and 
thereby cool the gas, hence restoring the temperature gradient. 
Sustained conduction and advection could further increase the 
central entropy. However, the timescale for such a process is 
long in comparison to radiative cooling timescales (< 1 Gyr) 
as well as cluster assembling and merging timescales which 
questions the possibility of this mechanism to boost the central 
entropies to values of K e ,Q ~ 600keVcm 2 wh ich are at the tail 
end of the distribution dCavagnolo et alj|2009l) . 

Because blazar heating does not inject large amounts of en- 
tropy at z > 2, and blazars do not efficiently heat high den- 
sity regions (1 +5 > 10), objects that have already collapsed by 
that time or shortly thereafter will not be significantly affected 
by blazars (barring the effect of clustering bias at early red- 
shifts). In our scenario, these early-forming groups evolve into 
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CC systems at the present epoch that potentially need to be sta- 
bilized by a self-regulated feedback process, e.g., provided by 
the radio mode of AGNs. However, the subset of groups that 
forms after z — 1 can be severely affected. If such a group is 
viewed shortly after its formation, it should still exhibit the ele- 
vated core entropies associated with the prior blazar heating. If 
such a late-forming group has a cooling timescale long in com- 
parison to the interval between cluster/group mergers, merger 
shocks can gravitationally reproce ss the entropy cores and am- 
plify th em b y a factor of up to five dBorgani et al.l2005l see also 
Section ing . These late-forming groups would then evolve into 
NCC systems^ For sufficiently late-forming clusters that ex- 
perience a series of fast successive merger events and, hence, 
avoid substantial cooling phases, we expect that gravitational 
reprocessing should boost the central core entropy in the most 
extreme case from ~ lOOkeVcm 2 to ~ 600keVcm 2 , allowing 
for reprocessing of the blazar-heated entropy floor due to grav- 
itational heating. Typically, however, blazar heated entropies 
at the time of turnaround for late-forming groups/clusters are 
~ 50keVcm 2 . Allowing for modest cooling periods in be- 
tween mergers should yield smaller median core entropy values 
of ~ (100- 150)keVcm 2 after accounting for gravitational re- 
processing. These esti mates compare favorably with observed 
values of NCC clusters dCavagnolo et al.ll2009l) . 

We point out that it is very natural for systems with a uni- 
modal distribution in core entropy values following group for- 
mation to evolve into a bimodal distribution today as conse- 
quence of both the cooling instability and the gravitational re- 
processing of (temporally increasing) elevated entropy cores. 
The observed CC and NCC cluster populations are centered 
upon the two attractor solutions that a galaxy cluster can evolve 
into. In the case of CC clusters, the core evolution is driven by 
the well-known overcooling problem, encountered in cosmo- 
logical simulations of galaxy clusters. Below a critical core 
entropy, purely hydrodynamical mergers are incapable of dis- 
rupting ac^mpac^CCsystems and transforming it into an NCC 
object dPoole et al.l2 008). For NCC clusters, the core evolution 
is driven by the rapid (in comparison to f coo i, which effectively 
sets the critical core entropy value within a given epoch) suc- 
cession of mergers, which has been demo nstrated to further ele - 
vate the core entropy values substantially dBorgani et al.l2005b . 
Interestingly, this solution is not a runaway solution; instead, 
the gravitational bootstrapping of blazar preheated entropy ad- 
justs to the system size. Hence, in this picture, the core entropy 
should never be able to exceed the entropy at the virial radius 
(according to virial arguments) and most likely reach only val- 
ues that are a fraction of that (at least for cluster systems) due 
to radiative cooling and the modest blazar preheated entropy 
values in comparison to early preheating models. 

To conclude, we demonstrate explicitly that the core entropy 
values in CC clusters have a very weak correlation with the me- 
chanical energy and power of X-ray cavities inflated by AGNs. 
Hence, PdV work done by these expanding cavities is trans- 
ferred inefficiently to the surrounding medium. This strongly 

15 Using Chandra observations, there has been a claim that incidence rate 
of cool core clu sters at redshifts z > 0-5 is much smaller than their fraction 
at low redshifts ( Vikhlinin et al. 20071). We note that the development of the 
characteristic cuspy X-ray brightness profiles for cool cores at z = requires 
the cooling time to be considerably shorter than the formation time — a crite- 
rion that is often not fulfill ed at the (high) r edshif t in question. Second, in order 
to construct their sample, Vikhlinin et al. 1 2007) discarded regions with AGN 
emission which can be correlated with cool cores, potentially biasing their ab- 
solute cool core rates low. Hence, this observation is not in contradiction with 
our proposed scenario. 



suggests that while AGN feedback seems to be critical in sta- 
bilizing CC systems, it cannot transform CC into NCC sys- 
tems (at least on the buoyancy timescale). With this evidence it 
seems even more pressing to pursue alternative solutions such 
as the presented blazar heating scenario in combination with 
gravitational reprocessing that provides a plausible scenario for 
the observed CC/NCC bimodality. Future cosmological hydro- 
dynamical simulations that include the effect of clustering bias 
of blazar heating are needed to study these considerations in 
greater detail 

2.4. Implications for the Sunyaev-Zel'dovich power spectrum 

The thermal SZ effect provides a direct probe of the gas 
properties of groups for z > 0.5. The SZ effect arises 
from CMB photons that inverse Compton scatter off ther- 
mal electrons within the hot plasma in galaxy clusters and 
groups, producing a localized perturbation to the CMB 
spect rum dSunyaev & ZeFdovichll 19721 : ISunyaev & Zeld ovichl 
119801) . The thermal SZ effect directly measures the thermal 
electron pressure in the gas and has the important property 
that its amplitude is independent of redshift. The pressure 
fluctuation spectrum of unresolved groups and clusters dom- 
inates the CMB power spectrum on angular scales smaller than 
3' (corresponding to a multipole moment I ~ 3000) and half 
of the SZ power spectrum signal at £ ~ 3 000 comes from 
groups withMsnn < 2 x 10 14 M Q and z > 0.5 dTrac et al.ll20l"il 
iBattaglia et aljl20TlT) . At these scales, the SZ power spectrum 
depends on the square of the Fourier transform of the average 
pressure profile of clusters/groups; a more concentrated pres- 
sure profile implies more power, a smoother one less. 

A population of groups/clusters with high core entropies, 
Kq = (50-100) keVcm 2 , implies a smoother pressure core dis- 
tribution than if the cores had cooled since formation and de- 
veloped a more concentrated entropy profile. This is very sim- 
ilar to the effect of AGN feedback which injects entropy into 
the core of groups, sm oothing out the resulting pressure profile 
dBattaglia et al.l 120101) . Adopting empirically motivated "uni- 
versal" pressure profiles constrained by X-ray observations, 
the peak amplitude of the SZ power spectrum is reduced for 
NCC pressure profiles as com pared to those of CC clusters 
(Efst athiou & Migliaccioll2"01 ll) . Hence, we expect the effect 
of blazar heating to result in a suppression of the SZ power 
spectrum for scales £ > 2000, which probe the pressure profile 
mostly inside R^oo- An abundant population of late-forming 
groups with high /To-values (that has not been taken into ac- 
count in any numerical mod eling o f the SZ power sp ectrum 
so far, e.g., by IBattaglia et al.ll2010l or lTrac et al.ll201 lb would 
furthermore reduce the thermal SZ power spectrum in compari- 
son to these numerical approaches. This has potentially impor- 
tant observational consequences since angular scales around 3' 
are the sweet spot for current telescopes measuring the high-£ 
CMB angular power spec t rum, e.g., the South Pole Telescope 
(SPT: [Lueker et al.l I20T0I: IShirokoff et al.l IMlt iKeisler et al.1 
1201 U IVanderlinde et al.l 120101) and the Atacama Cosmology 
Telescope (ACT; see, e.g.. iFowler et al.l 120101 ; iDunklev et all 

16 A detailed prediction of the core entropy distribution of clusters at a given 
mass (or temperature) depends on the merger history of all objects at that mass 
(to understand which number fraction of clusters was channeled into the cool- 
ing branch due to comparatively fast radiative cooling), the distribution of the 
departure times of the gas that ended up in the core from average densities (to 
address the exact magnitude of blazar heating for the population of clusters 
rather than individual systems), and quantifying the efficiency o f gravitational 
re-pro cessing for the entire population of groups and clusters ffiorgani et al.l 
12005) only simulated a group and a cluster for different variants of physics). 
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12011 iMarriage et al.ll20ll . 

In addition to the astrophysical dependence, the amplitude 
of the SZ power spectrum also depends very sensitively on 
cosmological parameters that are responsible for the growth 
of structure, Q oc al ' 9 {£lbh) 2 , where the rms amplitude of 
the (linear) density power spectrum on cluster-mass scales is 
denoted by a$. Hence, there is an interesting degeneracy in 
the amplitude (and shape) of the SZ power spectrum between 
the cosmological information (dominated by a&) and the as- 
trophysical information contained in the average pressure pro- 
file. Currently, numerical models of the S Z power spectrum are 
consistent with the d ata at the 1 a level (iDunkley et al.ll20l"lT : 
IShirokoff et al.ll20lll) . However, afte r allowing for a substan - 
tial signal from patchy reionization dlliev et alj 120071 [2008), 
which boosts the kinetic SZ effect and hence the total SZ sig- 
nal, the power predicted by these models becomes uncomfort- 
ably high. Suppression of power due to blazar heating rep- 
resents a promising mechanism by which to reconcile the ex- 
pected and observed SZ power spectrum. 

Unfortunately, the impact of blazar heating is degenerate 
with energy injection from AGN feedback in the dense cores 
of groups at early times (where there are also very little con- 
straints). Here, we sketch a promising idea on how to discrim- 
inate between the effects of AGN feedback and blazar heating 
on the SZ power spectrum. In particular, AGN feedback and 
blazar heating vary as a function of cluster mass and redshift. 
Their effects on the pressure profile of groups and clusters and, 
hence, their imprint on the SZ signal also vary as a function 
of cluster mass and redshift. By cross-correlating CMB maps 
(where foreground components and primary anisotropics have 
been subtracted) with deep optical redshift surveys that are sep- 
arately binned in redshift and cluster mass, i.e., optical richness 
estimators, we can perform SZ tomography. Su ch tomographic 
SZ po wer spectra (similar to Figures 7 and 8 of lBattaglia et alj 
l20Tll for the case of AGN feedback) will enable us to derive a 
redshift and mass-dependent mean cluster pressure profile. Po- 
tentially, this could disentangle the different non-gravitational 
energy injections of AGN feedback and blazar heating in clus- 
ters. However, large cosmological hydrodynamical simulations 
are required to obtain detailed predictions of either process 
which shall be subject to future work. 

3. STRUCTURE FORMATION AND DWARFS 

The entropy injected into the IGM by blazars not only mod- 
ifies the structures of groups and potentially clusters, but also 
has an observable effect upon cosmological structure forma- 
tion. Heating the IGM produces higher IGM pressures, which 
in turn suppress the gravitational instability on sufficiently 
small scales. Thus, there is a characteristic length scale, and 
hence characteristic mass (Mc), below which objects will not 
form. The particular value of this critical mass depends upon 
how structures form in practice, and generally requires a fully 
nonlinear study of structure formation (e.g., that provided by 
large-scale numerical simulations). Nevertheless, we may esti- 
mate the relevant characteristic length and mass scales via lin- 
ear perturbation theory. 

A very rough idea of the impact of blazar heating may be 
derived from the Jeans wavenumber, obtained by balancing the 
sound crossing and free-fall timescales: 



kjia) : 



c s {a) 



A-KGp(a). 



(9) 



factor, a = 1/(1 + z), p{a) = il,„(a)p cr {a) is the mean total mass 
density of the universe, c s (a) = y/5kT(a)/3pm p is the linear 
sound speed (in which T denotes the temperature at mean den- 
sity and p = 0.588 denotes the mean molecular weight for a 
fully ionized medium of primordial element abundance). In 
a static background universe, perturbations on scales smaller 
than lira/kj are stable, i.e., the ambient gas pressure is suffi- 
cient to counteract gravitational collapse. As a result, cooling, 
fragmentation, and star formation are suppressed within ob- 
jects of mass less than the Jeans mass, Mj, 
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The Jeans mass associated with a universe which has been 
heated by TeV blazars, Mjb\ azm -, is necessarily larger than one 
which has only undergone photoionization heating, My. p h to by 
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where we used a temperature ratio of 7bi azar /r p h to > {7, 14} 
at the present epoch for the standard and optimistic model, 
respectively (see Figure Q] and the discussion surrounding the 
standard and optimistic blazar heating histories). That is, 
blazar heating increases the mass of the smallest collapsed ob- 
jects by more than an order of magnitude with the exact value 
depending on the adopted numbers of blazars contributing to 
the heating. 

While the Jeans mass provides a way to estimate the impor- 
tance of blazar heating, Mj typically exceeds Mq by up to an 
order of magnitude because it neglects the Hubble expansion. 
That is, it fails to account for the time required for the pres- 
sure to influence the evolving gas distribution. A more rigorous 
treatment of linearized density perturbations in a baryon-dark 
matter fluid finds that the overdensity in dark matter, Sd(t,k), 
and the overdensity in the baryon distribution, 5b{t,k), both of 
which are functions of time and comoving wavenumber, are 
related by 

6 b (t,k) 



= 1 



k 2 

F 



5 d (t,k) 

for some kf dGnedin & Hui|[l998l) . The "filtering scale", as- 
sociated with kf, defines a size below which baryonic pertur- 
bations are smoothed despite the growth of background dark 
matter perturbations. The filtering wavenumber can be related 
to kj by 
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1 
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where D+(t) is the li near growth fu nction and is dependent 
upon the cosmology (Gnedirj [20001) . While kf is related to 
kj, at any time those can be very different since kf is an inte- 
gral over the past evolution of the Jeans scale, weighted by the 
appropriately scaled growth function. 

The associated mass scale, defined in analogy with the Jeans 
mass, is 

„ < s 47r -, s ( 2 ™ X 

M F (a) = — p(a) — — 
3 \k F {a) 



(14) 



where for convenience, we have introduced the cosmic scale 



(the details of how this is computed in practice are collected 
in Appendices lAl andlBl . To understand the physics underlying 



THE COSMOLOGICAL IMPACT OF BLAZAR TEV EMISSION III 



15 



10 l; 



© 10'"- 



1 + 8= 1,2-™= 10 



10 s 



10" 




a 10 r 

3 - 

sf it 



linear theory 
non-linear theory 
optimistic blazar h. 
standard blazar h. 
only photoheating 



10 



10' 



10 ,u - 



10 s 



10° 



I 10 



1+8 = 0.5 >W =10 




linear theory 
non-linear theory 
optimistic blazar h. 
standard blazar h. 
only photoheating 



10 



1+z 



1+z 



Figure 7. Redshift evolution of the filtering mass, Mf, for the cosmic mean density, <5 = (left) and for a void with mean overdensity, 8 = —0.5, (right). We contrast 
Mf in the standard cosmology that employs only photoheating (solid) to the case of blazar heating in our standard model (dashed) and optimistic model (dash-dotted). 
In the bottom panels, we show the ratio of Mf in our respective blazar heating models to those without blazars. To estimate the effect of nonlinear structure formation 
on the filtering mass, we compare th e linear theory Mf (blue) to the nonlinear theory Mf (red) where we used a correction function derived from hydrodynamic 
simulations by Okamoto et al. 1 2008). The dotted line represents the analytical solution for Mf which is only valid before reionization for a power-law temperature 
evolution. 



the filtering mass, it is instructive to connect the filtering mass 
to the entropy of the IGM. To this end, we use the definition 
for the entropy K of Equation (fTJ. As we show in AppendixlAl 
the filtering scale Xf = lua/kf is an integral over the entropy 
evolution, 

'^Af^^r^, (15) 

4(a) D + (a)J 'a«E(a')J a , a"iE(a»y 

. 5 ( 3 a,, \ 2/3 

Ao = 3{^GH ) ■ (16) 

Here E(a) = H(a)/H{) is the dimensionless Hubble function. 
Since Mf oc IcJ 3 , the linear filtering mass is also determined by 
the entropy evolution at mean density and appropriately scaled 
with the linear growth function, Hubble function, and cosmic 
scale factor. Since the integrand in Equation ( fl5l ) is a positive 
definite function, the linear filtering mass increases monotoni- 
cally with time. 

We expect the strong increase of the entropy due to blazar 
heating to drive a comparably strong increase in Mf, though 
slightly delayed due to the effective weighting function in M F . 
Physically, this implies that the entropy floor delivered by 
photo- or blazar heating at densities around the cosmic mean 
is then conserved by adiabatic compression during structure 
formation and directly translates to a linear filtering mass Mf. 
However, during this nonlinear process of structure formation, 
the entropy can be additionally augmented through dissipation 
of gravitational energy in cosmic structure-formation shocks or 
decreased through radiative cooling. The competition of these 
two processes determines which will dominate and in turn sets 
the critical mass Mq which is required for the condensation of 
gas into a dark matter halo and eventually the formation of a 
galaxy. 

As with Mj, Mf typically exceeds the Mq observed in sim- 
ulations by a substantial factor. Whether a halo can accrete 
gas is determined by the gas temperature (or equivalently en- 
tropy) at the virial radius which can only be obtained through 
hydrodynamic cosmological simulations. Estimated values for 



this critical mass using cosmological simulations of nonlinear 
structure formati on with only photoheating yield Mg(z = 0) = 
6.5 x IOV'Mq dHoeft et al.l2006t lOkamoto et al.l2008l) — one 
order of magnitude sm aller than the linear analog Mf at z = 
( lOkamoto etal.ll2008l see also Figure|7]). This discrepancy be- 
comes less for increasing redshift leading to reasonably good 
agreement at high redshift (z > 6) between Mf and Mq- To 
model the nonlinear behavior of Mf, we can introduce a correc- 
tion factor C(z) = Mc(z)/M F (z) = (l + z) u /H-8, where M F (z) 
has been modeled from the temperature evolution at mean 
density in the cosmological simulations and includes weakly 
nonlinear aspects o f structure formation (lOkamoto et alj |2008; 
iMaccio et al.ll2010T) . However, validating this approximation 
for Mq will ultimately require numerical simulations that in- 
clude blazar heating. 

Figure [7] shows Mf as a function of redshift for a void and 
the cosmic mean when various heating mechanisms are consid- 
ered (which a grees for the standard cosmology with Figure 5 of 
Gnedin 2000). Note that for all of the thermal histories we con- 
sidered, the resulting Mf nicely follows the expected analyti- 
cal form over its full range of validity prior to recombination 
(which necessarily implies an Einstein-de Sitter universe, see 
the derivation in Appendix|B]). While blazar heating increases 
Mf in voids as early as z ~ 3, perturbations at the cosmic 
mean are affected slightly later (z ~ 2.5). The extra suppres- 
sion of the formation of small galaxies due to blazar heating 
amounts to a mass suppression factor of ^25-70 for voids and 
~ 6.5 — 15 for the cosmic mean today. The range indicates the 
uncertainty of the blazar heating models and the larger value 
belongs to the model that matches the inv erted temperature - 
density relation in the Lya forest found by IViel et all ( 120091) . 
i.e., the optimistic model. The resulting estimates for the non- 
linear characteristic mass Mq are also shown in Figure [7] with 
red lines for our models with and without blazar heating. Ap- 
parently the increase of Mq due to blazar heating is able to 
counteract the suppression of Mf in nonlinear theory. We note 
that our estimate for Mq is slightly decreasing for z < 1 in the 
model that employs only photoheating. This decrease is an ar- 
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tifact of our linear theory Mf which by construction does not 
model nonlinear aspects of structure formation such as shell 
crossing or formation shocks that raise the temperature of some 
patches already at mean density. 

We now turn to two outstanding problems in cosmological 
structure formation that the recent blaz ar he ating may help 
to addre ss: t he missing satellite (Section 13. U and void dwarf 
(Section [i.2\ problems. Finally in Section [331 we discuss how 
blazar heating may naturally address some of the difficulties 
with WDM cosmologies of galaxy formation. 

3.1. Dwarf satellites in the Milky Way 

The heating due to blazars provides an additional mecha- 
nism to suppress the formation of dwarfs. Unlike photoion- 
ization models, which typically invoke the heating at reioniza- 
tion, blazar heating provides a well defined, time-dependent 
suppression mechanism, with the suppression rising dramati- 
cally after z ~ 2. This can be seen explicitly by the steep in- 
crease in filtering mass for these redshifts in Figure [7] In addi- 
tion, due to the homogeneous nature of the heating, i.e., a con- 
stant volumetric rate that is independent of density, the heating 
from blazars suppresses structure formation most efficiently in 
the low-density regions responsible for late-forming dwarf ha- 
los due to their negative bias, where the energy deposited per 
baryon is larger. 

In fact, the star formation histories of dwarf galaxies provide 
a strong constraint upon the magnitude of blazar heating: if 
dwarf galaxy formation is suppressed by high-energy gamma- 
ray emission from blazars, all dwarf star formation histories 
must begin prior to z — 2, roughly the redshift at which blazar 
heating becomes significant. Over the past decade most of Lo- 
cal Group dwarfs have been observed with the Hubble Space 
Telescope shifting the focus from counting dwarfs to resolving 
the individual stellar population within these objects. Thus, it 
has become possible to construct detailed color-magnitude di- 
agrams of every dw arf galaxy and therefore their detailed star 
formation historie s (Dolphin et al. 2005c iHoltzman et al.l l2006; 
lOrban et "ail 120081) . While the data shows a great variety of 
star formation histories — some continuous, some bursty, some 
truncated — they all have in common that they extend beyond a 
lookback time of lOGyr (corresponding to z = 2). There is no 
known case of a dwarf galaxy that formed its first set of stars 
after z — 2. The fact that there exists an old stellar population 
in every dwarf in the Local Group provides a very important 
test that our model of blazar heating successfully passes. 

The degree to which blazar heating can suppress the num- 
ber of dwarf galaxies depends upon the redshifts at which 
they are typically formed. By combining high-resolution N- 
body simulations of the evolution of Gal axy-sized halo with 
semian alyti c models of galaxy fo rmation, IMaccio & Fontanot] 
(120101) and IMaccio et all (120101) have inferred the statistics, 
formation timei' l and accretion histories of the Milky Way 
satellites. Their findi ngs are in good agreement with recent 

^2009tlMunoz et al.ll2009b 



work by other groups (iKoposov et a 



IGuo et al.ll2QTol: iBusha et alj|2010t iFont et alj|201 lb and thus 
can be regarded as representative for this model class. They 
find that the distribution of formation times is bimodal as a re- 
sult of the suppression of hot gas accretion in low-mass halos 

17 IMacc io et al. (20lT|) defined the formation time of dwarfs as the redshift 
where the progenitors of dwarfs that end up within the Milky Way halo today 
exceeded a virial temperature of Z200 = 10 4 K such that H I cooling of the gas 
becomes possible in these halos. That virial temperature corresponds to a mass 
threshold of M 2 oo ~ 10 9 M Q at z = 1 and 10 8 M Q at z = 10. 



due to the photoionization background; while 2/3 of today's 
satellites form at redshifts ranging from 3 < z < 12, 1/3 of all 
satellites form late at z < 3 with most of them at z < 1.5. Be- 
cause their formation is marked by the first time H I cooling can 
effectively form the first stars, their stellar populations are nec- 
essarily younger than their formation time of z < 1 .5 in direct 
conflict with observed ages of the stellar populations in all of 
the Local Group dw arffl ( D olphin et alj2005tlHoltzman et al 



2006; lOrban et alj 120081) . The model of IMaccio & Fontanot 



2010ft that successfully reproduced the satellite luminosity 
functi on of the Milk y Way assumed the linear theory filtering 
mass (lGnedinll200d see solid blue line in Figure IT) which was 
shown to significantly overproduce the characteristic halo mass 
scale below which b aryons cannot condense and form stars 
( Ok amoto et alj|2008l see solid red line in Figure |7J. Adopting 
their less efficient nonlinear filtering mass formalism, the num- 
ber of satellite galaxies with intermediate magnitudes My ~ 10 
increases and creates a bump in the luminosity function which 
is clearly inconsistent with the data (Maccio et al.!2010l). 

Inspection of nonlinear estimates for the filtering mass Mf 
in Figure [7] shows that heating from blazars could prevent the 
condensation of baryons in halos of masses 10 10 M o at red- 
shifts z ~ 2 and likely up to 10 n M Q at z = 0. At mean den- 
sity, these nonlinea r estimates for Mf correspond to the origi- 
nal linear values of lGnedinl (120001) which have been shown to 
yield the observed dwarf satellite abundances a nd luminosity 
functions ( ISomervillell2002t IMaccio et al.ll2010l) . At the same 
time, the strongly rising Mf in the blazar heating models after 
z ~ 2 is able to suppress th e population of late-forming dwarfs 
(Ma ccio & Fontanotl 120101) that are in conflict with measured 
star formation histories of dwarf galaxies. Thus, not only can 
blazar heating potentially play a significant role in explaining 
the observed abundances of dwarf galaxies but also their old 
star formation histories. 

Finally, we argue that the redshift evolution of the blazar 
heating rate should manifest itself as stochasticity in the satel- 
lite luminosity function at fixed host halo ma ss. Host halos 
show a distribution of formation times (e.g., IWechsler et alj 
2002). This distribution is inherited by the host's satellite 
dwarfs which on average form earlier than the host. As demon- 
strated in Figure[T| blazar heating implies an entropy floor that 
dramatically increases after z ~ 2. Hence, the distribution of 
formation times of satellite dwarfs at fixed host halo mass re- 
sults in a distribution of pre-collapsed entropy of these dwarf 
halos. As a result, this leaves early-forming dwarfs relatively 
unchanged but suppresses the baryon fraction or even the for- 
mation of late-forming dwarfs. This results in different cool- 
ing histories for different dwarfs, which, in turn, might modify 
the stellar content that condenses out in these systems. Hence, 
blazar heating provides an apparently substantial, physically 
motivated scatter in the satellite luminosity function, or dis- 
tribution of mass-to-light ratio at fixed host halo mass. This 
physically motivated stochasticity has important implications 
for abundance matching techniques that need to be taken into 
account and complicates their use. We note that some of this 
stochasticity is already seen in mildly nonlinear theory (which 
captures some aspects of the formation history) and manifests 
itself as a scatter in entropy at large l+S in FigurefJJ 



18 This assumes star formation in primordial gas. However, other processes 
such as metal pollution from adjacent galaxies might cause stars to form earlier 
at lower virial temperatures of T200 < 10 4 K. 
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3.2. The Void Phenomenon and the Faint-End Slope of the 
Galaxy Luminosity and H l-Mass Function 

As outlined in Section [TT2l the "void phenomenon" is closely 
related to the substructure problem. Both show a strong dis- 
crepancy in the abundance of dark matter (sub-)halos and 
paucity of luminous dwarf galaxies that are thought to be 
hosted by these halos. Before we discuss in detail how blazar 
heating impacts on dwarf galaxies in voids, we review the latest 
discussion on the status of the problem itself which has recently 
been disputed to even exist. 

Using a halo occup ation distribution model approach, 
iTinker & Conrovl ( 120091) claim to have explained the problem 
as they find agreement between luminosity functions, nearest 
neighbor statistics, and void probability function of faint galax- 
ies. However, a very high resolution simulation of the local 
volume of 8 Mpc around the Milky Way predicts a factor of 
10 more dwarf halos than observed dwarf galaxies in mini- 
voids with sizes ranging from 1 to 4.5 Mpc, hence reinforcing 
the "void phenomenon" (iTikhonov & Klvpinll2009t) . While the 
agreement between theory and observations is good for dwarfs 
with masses M%x)m si 1O 1O M (corresponding to maximum cir- 
cular velocities of v c > 40 km s _1 )Fjit fails below, suggesting 
that lTinker & Conrovl (|2009) did not sample the relevant mass 
scales. Moreover, their analysis only demonstrated the self- 
consistency of the halo occupation distribution model and did 
not make a comparison between observations and the predic- 
tions of ACDM. 

The discrepancy of the luminosity function on dwarf scales is 
confirmed by other recent studies that compare the circular ve- 
locity function of H I observa tions with those obt a ined through 
dissipationless simulations dZavala et al.l 120091; IZwaan et al.l 
120101: iTruiillo-Gomez et al.ll20llb . At v c < 80 km s" 1 (corre- 
sponding to M2oom < 10 u Mq), the latest study finds a slight 
deviation between theory and observations which amounts to 
a significant overprediction of more than 10 times the number 
of observed systems at v c = 40 km s -1 . While completeness of 
the observed samples could be an explanation of some of the 
differences, it is unlikely that it can account for all of the ob- 
served effects since the blind H I sample of the HIPASS survey 
is thought to be co mplete down to M» , < 5.5 x 1O 7 M out to a 
distance of 5 Mpc (IZwaan et al.ll2010» . Since gas-rich galaxies 
dominate at the low mass end of the luminosity function, their 
sample should give an accurate measurement of the abundance 
of dwarfs if these galaxies contain enough neutral gas to be 
detected. 

While blazar heating significantly increases the entropy, and 
hence the associated filtering mass, at mean density, it does 
so more dramatically within the voids. To address how this 
substantially increased heating manifests itself in the num- 
ber of void dwarfs, it is instructive to compare the formation 
timescales of dwarf halos as a function of environment. Void 
dwarf galaxies are negatively biased and form later than their 
field or cluster analogs (halos of a given mass tend to be older in 
clusters and younger in voids for masses smaller than the typi- 
cal mass scale that is presently entering the nonlinear regime 
of perturbation growth). This is because for galaxies form- 

19 We quote the values corresponding to the lower bounds in their models 
that assumed a normalization of the matter power spectrum of og = 0.9. This 
seems to be a conservative choice when comparing to the critical mass thresh- 
old of M200111 > 6 X 10 s M© or equivalently v c > 35 km s _I which assumed a 
0% = 0.75. For consistency reasons with the literature, in Section [3] this sec- 
tion, we define the virial mass, M200M, as the mass of a sphere enclosing a 
mean density that is 200 times the mean density of the universe. 



ing on a large-scale underdense mode, more time has to elapse 
before these galaxies acquire enough overdensity to decouple 
from the Hubble expansion and subsequently collapse. Thus, 
the median redshift of formation for a dwarf galaxy with halo 
mass 2.4 x lO^M^ is Zf orm =2.1 within cluster and 1 .6 in voids 
(lHahn et alJl2007aT) . Moreover, the distribution of formation 
times of void galaxies is more sharply confined around the me- 
dian value than those for clusters, which exhibit a long tail o f 
formation redshifts, extending to z ~ 6 ( lHahn et al.l I2007b1) . 
Our estimate of Mq in underdense regions with 1 + S = 0.5 
demonstrates that blazar heating prevents the condensation of 
baryons in halos of masses of 1O 1O M at redshifts z < 2.4 and 
2.8 x 1O H M at z = (see Figure |7J. In particular, halos of 
2 x 1O 1O M can be suppressed after z = 1.8 which is earlier 
than the median formation redshift of these galaxies. Depend- 
ing on the exact definition of voids, this suggests that more than 
half of the galaxies with masses M < 2 x 1O 1O M , correspond- 
ing to a maximum circular velocity of v c < 45 kms -1 , can be 
suppressed or severely affected by blazar heating. 

Such a preheated entropy floor not only suppresses dwarf 
formation at late times (z < 2), but may also modify galaxy 
formation at the low-mass end. Assuming that low-mass halos 
are embedded in a preheated medium with an entropy floor of 
lOkeVcm 2 at z < 2 simultaneously matches data at the faint- 
end slope of the ga laxy luminosity function as well as of the 
H i-mass function dMo et alJ 120051) . This heuristic assump- 
tion almost exactly coincides with the predictions of our blazar 
heating models (see Figure [TJ. As a result of such a preheat- 
ing, only a fraction of the gas in a proto-galaxy region would 
be able to cool and be accreted into the final galaxy halo by the 
present time. If the accreted gas resides in the diffuse phase, 
it does not lose angular momentum to the dark matter, thereby 
possibly conti nuing to form larg e galaxy discs in low-density 
environments dMo & Madl2002h . 

In summary, the entropy floor and filtering mass due to blazar 
heating are dramatically increased in voids as a result of the 
constant volumetric heating rate, which leads to an inverted 
temperature-density relation in low-density regions (Paper II). 
In combination with the later formation epoch of dwarfs at 
these low densities, this implies a very efficient mechanism for 
suppressing void dwarf formation in collapsed dark matter ha- 
los. Hence, our model provides an elegant physica l solution to 
the void phenomenon described by Peebles! (120011) . 

3.3. Suppression of Dwarfs in Warm Dark Matter 
Cosmologies 

Recent dissipationless ACDM simulations produce not only 
far too many dark matter satellites but also show that the most 
massive subhalos in simulations of the Milky Way may be too 
dense to host any of its observed bright satellites with luminosi- 
ties Ly > 10 5 L Q , i.e., these massive satellites in simulations at- 
tain their maximum circular velocity at t oo small radii in com- 
parison to the observed dwarf satellites dBoylan-Kolchin et alJ 
1201 lblah . These dark subhalos have circular velocities at in- 
fall of 30-70km s" 1 and infall masses of (0.2-4) x 1O 1O M . 
In principle, this puzzle can be solved (or partially solved) in 
the following ways: by increasing the stochasticity of galaxy 
formation on these scales, by reducing the central (dark mat- 
ter) densities by means of very efficient and violent baryonic 
feedback processes acting on timescales much faster than the 
free-fall time, by assuming a total mass of the Milky Way 
at the lower end of the allowed uncertainty interval (i.e., ~ 
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8 x lO n M ) in combination with a shallower subhalo den- 
sity profile of the Ein asto form as measured in simulations 
dVera-Ciro et al.l 120121) . or by allowing these subhalos to ini- 
tially form with lower concentrations as would be the case, for 
example, if th e dark matter wer e made of warm, rather than 
cold particles dLovell et al.ll2012l) . 

In the limit of heavy WDM particles, e.g., the (sterile) neu- 
trino with m u c 2 > lOkeV that is created in the early uni- 
verse through mixing with an active neutrino, structure for- 
mation proceeds almost indi stinguishably from CDM for all 
current observational probes dSehak et alj|2006l) . For smaller 
masses, however, the free streaming of neutrinos erases all 
fluctuations on scales smaller than the free streaming length, 
which is roughly proportional to their temperature and in- 
versely proportional to their mass. Using the Lya forest 
power spectrum measured by the Sloan Digital Sky Survey and 
high-resolution spectroscopy observations in combination with 
CMB and galaxy clustering constraints still allows for a neu- 
trino with mass m v c 1 > 2.5 ke V (95% c.l.) that d ecoupled early 
while in thermal equilibrium dSeliak et alJ^OO^FI 

Dissipationless high-resolution simulations of the evolution 
of a Galaxy-sized halo have shown that (sub-)halos form and 
are accreted later onto the main halo in WDM scenarios com- 
pared to the standard CD M paradigm, due to the lac k of power 
on small scales in WDM dMaccio & Fontanotl20Tol) . In WDM 
scenarios with relatively low values of the particle masses, 
ra„c 2 = (2-5) keV, there are almost no halos with Zf orm > 11, 
suggesting that the fraction of late-forming dwarfs Zf olm < 1.5 
is increased over the CDM scenario. This late formation epoch 
of dwarfs reinforces the star formation history problem of Lo- 
cal Group dwarfs: these late-forming dwarfs contain young 
stellar populations, in direct conflict with the old ages of the 
stellar populations (r > 10 Gyr) in all of the Local Group 
dwarf s dDolphin et atll2005l: iHoltzman et"aDl2006l: IOrban etal.1 
120081) . 

However, these are precisely the objects that blazar heating 
most strongly affects, as the strongly rising filtering mass (or 
entropy floor) after z ~ 2 is able to suppress the population of 
these late-forming dwarfs. Hence, while in this scenario the 
free streaming of WDM erases power on scales smaller than 
dwarfs, blazar heating reconciles the theoretically expected and 
the observed star formation histories and alleviates standard 
objections to galaxy formation in WDM cosmologies. 

3.4. Impact on the formation ofL* galaxies 

In the previous sections, we argued that a blazar-heated en- 
tropy floor is able to suppress late dwarf formation in voids 
and the Milky Way as well as modify the thermodynamical 
profiles of galaxy groups and clusters. Hence, it is natural to 
ask whether there would be any effect of blazar heating on L* 
galaxy formation which represents the mass scale in between 
these two extremes. Blazar heating is not powerful enough 
to raise the mean temperature at 5 = of even the most ex- 
treme patches above 10 5 K (see Figure 9 in Paper II). The clas- 
sical criteria of galaxy formation are a short cooling time com- 
pared to the dynamical time and to the age of the universe, 
/cool < H~ l and t zoo[ < f dyn dRees & Ostrikerill977t Isiiklll977l: 
I Whit e & Rees l 1978b . which are easily fulfilled even in the pres- 

20 Since blazar heating dramatically changes the thermal history of the IGM 
(Paper II) and may be responsible for the inverte d temperature-density relation 
at z = 2— 3 inferred by high- redshift Lya studies (Bolton et al. 2008; Viel et al. 
2009; Puchwein et al. 201 1 ), it is not clear whether these limits on sterile neu- 
trino properties are weakened in the presence of blazar heating. 



ence of blazar heating (see, e.g., Figure 1 in iRees & Ostrikefl 
11977b . This implies that these galaxies can radiate away the 
additional entropy that the gas attained prior to collapse due 
to blazar heating within a free-fall time. It is, however, inter- 
esting to speculate whether the blazar-heated high-entropy gas 
at low redshift z < 1 has any impact on the late accretion of 
gas into the hot reservoir of baryons from which gas cools and 
fuels the late time star formation. Higher entropy gas should 
shock further out in the halo than pre-cooled gas of lower en- 
tropy which is denser and can provide a larger ram pressure. 
Blazar heated gas has high entropy and is more dilute. Hence, 
it is more easily torqued (by dissipative processes or magnetic 
fields) and therefore change its angular momentum distribu- 
tion. Thus, it is implausible that blazar heating dramatically 
changes the ordinary mode of galaxy formation, but we might 
anticipate blazar heating to starve the late-time accretion and 
potentially slow down subsequent star formation. These ideas 
are subject to verification by numerical simulations of galaxy 
formation. 

4. CONCLUSIONS 

TeV blazar heating results in a dramatic increase in the en- 
tropy of the IGM following He II reionization around z ~ 3.5. 
Since the IGM entropy evolution is critical for the formation 
and structure of collapsed objects, blazars heating has a signif- 
icant impact upon both. We have identified two mass ranges 
(or classes of objects) for which the TeV blazar-induced en- 
tropy floor should have a substantial effect. Galaxy groups and 
clusters, which are forming near the peak entropy injection rate 
(z ~ 1) and exhibit core entropies that are comparable to that 
implied by blazar heating; and dwarf galaxies, which are sus- 
ceptible to the rapidly rising entropy floor generated by blazars. 
Below, we describe the consequences for each in more detail. 

Galaxy groups and clusters. Immediately after formation, 
groups at fixed mass should have a continuous distribution of 
core entropy values, depending on the formation redshift and 
the temporally variable heating mechanism. The fate of these 
groups is determined by the ratio of the cooling time, f coo i, to 
the timescale between cluster/group mergers, f me rger- If this ra- 
tio is smaller than unity, the group can radiate the elevated core 
entropy away and evolve into a CC which survives the succes- 
sive hierarchical growth. Alternatively, if f coo i > f mer ger, merger 
shocks can gravitationally reprocess the entropy cores and am- 
plify them. Those groups can then evolve into NCC systems. 
Hence, it is not necessary to produce all of the observed cen- 
tral entropy in the IGM before collapse, but it is also possible 
to achieve this through gravitational heating, provided there is 
a certain minimum entropy delivered by some putative heat- 
ing process corresponding to a minimum f coo i. An increasing 
entropy floor also implies an increasing cooling time, hence 
the cluster-averaged f coo i /fm erger increases. It follows that sys- 
tems that evolve into CC systems today are on average early- 
forming, i.e., old systems. In contrast, NCCs are on average 
young systems. 

We argue that systems with a unimodal distribution in core 
entropy values after group formation should naturally evolve 
into a bimodal distribution. The reason for this are the two 
attractor solutions of the group/cluster system, driven by the 
cooling instability and gravitational reprocessing of (tempo- 
rally increasing) elevated entropy cores, resulting in the ob- 
served populations of CC and NCC systems, respectively. Such 
an elevated entropy core level in groups might explain a pop- 
ulation of X-ray dim groups with low gas fractions. We show 
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that the core entropy values in CC clusters have a very weak 
correlation with the mechanical energy and power of X-ray 
cavities inflated by AGNs. Apparently, PdV work done by the 
expanding cavities is an inefficient heating process that does 
not generate much entropy. This strongly suggests that while 
AGN feedback seems to be critical in stabilizing CC systems, 
it cannot transform CC into NCC systems (at least on the buoy- 
ancy timescale). 

Our blazar-induced entropy history seems to be well matched 
for the formation times of today's groups, but at first sight less 
so for clusters which form in highly biased regions through 
mergers of groups. Those had to form even earlier when the 
entropy for the average IGM was still rising with typical val- 
ues at z — 1 of Kq ~ (25-50)keV cm 2 . However, two effects 
positively interfere to counteract the apparently smaller effect 
of blazar heating in clusters. First, the mass accretion rate is 
larger for larger systems and at higher redshifts. This suggests 
that the earlier forming group progenitors of clusters can toler- 
ate a smaller blazar-heated entropy floor after collapse which 
will then be gravitationally processed at a faster rate and able 
to counteract the smaller cooling timescales. Second, blazars 
also turn on first in highly biased regions, and thus the IGM 
in the vicinity of clusters should experience blazar heating ear- 
lier than low-density regions. We speculate that this effect in 
combination with the faster gravitational reprocessing of an el- 
evated entropy core in denser regions would help in propagat- 
ing the effect of blazar heating to the scale of massive clusters. 

Changing the thermodynamic structure of groups and po- 
tentially clusters has also an impact on the SZ power spec- 
trum which is sensitive to cosmological parameters such as erg 
and the thermal pressure profile. An increased core entropy 
level implies a smoother pressure profile and hence decreases 
the power in the SZ power spectrum on scales smaller than 
3' which are probed by current experiments such as SPT and 
ACT. Lowering the astrophysical signal can allow for larger 
values of erg, potentially reducing some (minor) tension with 
the data especially after allowing for a contribution due to 
patchy reionization. 

Dwarf galaxy formation. We demonstrate that the redshift- 
dependent entropy floor increases the characteristic halo mass, 
Mc, below which dwarf galaxies cannot form by a factor of 
6.5-15 for the cosmic mean (6 = 0) and by a factor of 25-70 
for voids (S = -0.5). The range indicates the uncertainty in the 
number of blazars that contribute to the heating and the upper 
envelope matches the observations of an inverted temperature- 
density relation in the Lya forest. The increase of Mc pre- 
vents the formation of late-forming dwarf galaxies (z < 2) with 
masses ranging from 10 10 to 10 11 M Q for redshifts z ~ 2 to 0, 
respectively. This may resolve the "missing satellites problem" 
in the Milky Way, i.e., the low observed abundance of dwarf 
satellites compared to CDM simulations. It also brings the 
observed early star formation histories of Local Group dwarfs 
into agreement with galaxy formation models that predicted a 
population of late-forming objects, in conflict with the data. 
At the same time, it provides a plausible explanation for the 
"void phenomenon" which is the apparent discrepancy of the 
number of dwarfs in low-density regions in CDM simulations 
and the paucity of those in observations. Blazar heating sup- 
presses the formation of galaxies within existing dwarf halos 
of masses < 3 x 10 10 M Q with a maximum circular velocity 
< 60 km s -1 for z < 2. Additionally, the phenomenology of 
such a preheating mechanism matches heuristic assumptions 
that were adopted to match the faint-end slope of the galaxy 



luminosity function as well as of the H 1-mass function, in par- 
ticular for low-density environments. 

We conclude that the presented scenario of blazar heating 
holds the promise for solving some of the most outstanding 
problems in high-energy gamma-ray astrophysics, the IGM as 
probed by the high-redshift Lya forest, the formation of galax- 
ies, and clusters of galaxies. At the same time, it provides an 
astrophysical solution to these problems such as the "missing 
satellite problem" or the "void phenomenon" which have been 
claimed to require new physics beyond the concordance cos- 
mological ACDM model. 
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APPENDIX 

A. COMPUTING THE FILTERING MASS GENERALLY 

Here we collect the relevant details for computing the filter- 
ing mass, Mf, associated with a particular cosmological and 
thermodynamic evolution of the IGM. Recall that the filtering 
wavenumber is related to the Jeans wavenumber via 

1 - 1 r^ flV) ^+w)z3 +( o 



m) d + (o 



kj(t') 



dt' 



/,.w m 

We can simplify Equation ( TT3T > by noting that the linear growth 
function obeys the following equation 



D + (t) + 2H(t)D + (t) = 4irGpD + (t), 
in which the Hubble function is given by 
H\a) 



Hi 



■ E z (a) = a~ i fl m + a\\ - Q. m - VL a ) + VL 



A ■ 



(Al) 



(A2) 



Substituting the integration variable in Equation ( U~3l l by the 
scale factor a, we obtain 
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Following ICarroll et al.l dl992l) . the linear growth function can 
be computed by 
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For any given temperature evolution of the IGM (which enters 
via the sound speed), we can co mpu te a redshift evolution of 
the filtering scale from Equation dA3l >. It is convenient to define 
a filtering mass by analogy with the Jeans mass: 



4tt _ ( 2na 
M F (a) = —p(a) — — 
3 \k F (a) 



(A5) 



Finally, using a definition for the entropy of Equation ([T), we 
can derive the filtering scale X F = 2-Ka/k F in terms of the en- 
tropy, 

k 2 F (a) D + (a)J 'a ,3 E(a')J a , a" i E(a"Y y 



3 \SttGHoJ 
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B. FILTERING MASS IN AN EINSTEIN-DE SITTER UNIVERSE 

For sufficiently early times, all matter-dominated 
Friedmann-Lemaitre model universes asymptotically ap- 
proach an Einstein-de Sitter universe (z > 2 for our ACDM 
universe), for which the cosmological constant A = and 
the curvature = 0. For such a universe, the growth func- 
tion D+( a) oc a which allows us to considerably simplify 
Equation iA3i , yielding 
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Note that we replaced the unphysical lower integration limit 
by the recombination scale factor, a flmin , as baryon perturba- 
tions can only start to g row after recombination. Employing 
the relation of Equation ( IA5t we can rewrite this to give 



M 2 J\a)= - I da'M 2 J\a') 



a 



(B2) 



After recombination, the residual electron density couples 
the gas temperature still to that of the CMB via Compton in- 
teractions. Hence, we expect the gas temperature to scale as 
T oc a~ l . At z — 150, the Compton interaction rate drops below 
the Hubble expansion rate such that the gas experiences adia- 
batic expansion with T oc a~ 2 . Hence for the time after matter- 
radiation equality, where p = po a -3 , we can write in general 



T(a) 
To 
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This definition for the Jeans scale implies a Jeans mass at a 
fiducial scale factor, ao, of 



Mj.o 
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If we substitute Equation ( IB3b into Equation dB2l ). we can ob- 
tain an analytical solution for the filtering mass before reion- 
ization, 



M F (a) = Mj X )- 



da 1 



1/2' 
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Evaluating this integral for the two regimes before and after 
the freezeout of Compton interactions of the gas with the CMB 
photons, yields the following analytic solutions, 
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